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Preface 


Developments  in  high-energy  laser  candidates  are  occurring  at 
a rapid  pace.  As  the  number  of  candidates  for  the  Department  of 
Defense  High-Energy  Laser  Program  increases,  an  evaluation  or  each 
candidate  in  detail  is  more  difficult.  Therefore,  the  need  for  simple 
analyses  of  the  candidates  exists.  These  evaluations  need  to  give 
approximate  results  for  each  laser  so  that  research  can  he  directed 
toward  the  more  promising  prospects. 

This  paper  is  an  attempt  to  make  a simple  analysis  of  the  HgCl 
laser  which  is  a new  candidate  for  the  Department  of  Defense  High- 
F.nergv  Laser  Program.  Several  of  the  reaction-rate  constants  for  the 
reactions  in  the  HgCl  laser  are  computed  by  using  simple  theories 
applicable  to  these  reactions.  Some  reaction-rate  constants  are 
computed  by  using  experimental  data.  The  model  that  is  constructed  from 
the  reactions  and  reaction-rate  constants  is  analyzed  by  both  an 
analytical  solution  and  a computer  program  solution. 

Several  people  have  provided  me  with  assistance  and  aid  throughout 
this  study.  Without  the  contact  I have  had  with  these  people,  this 
study  would  lack  some  key  elements.  In  particular,  I want  to  thank 
Dr.  M.  R.  Flannery  of  the  Joint  Institute  for  Laboratory  Astrophysics, 
Boulder,  Colorado,  for  sending  me  copies  of  his  calculated  data  on 
Hg+  and  Cl"  ionic  recombination.  I want  to  thank  Pr.  L.  Schlie  of 
the  Air  Force  Weapons  Laboratory  for  several  interesting  conversations. 

A great  debt  of  gratitude  is  owed  to  A.  M.  Hunter  of  the  AFIT  Physics 
Department,  my  thesis  advisor,  for  the  invaluable  help  he  gave  me 
throughout  my  research.  I want  to  thank  Ms.  Sharon  Flores  for  typing 
the  final  copy  of  this  thesis.  Finally,  a special  note  of  thanks  to 
my  wife,  Kathy,  for  her  patience  and  many  hours  spent  typing  the  rough 
draft  of  this  thesis. 
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REACTION-RATE  CONSTANTS  AND  SIMPLE  RATE  THEORIES 
APPLIED  TO  THE  MERCl  RY-CHl.ORIDl  LASER 


I • Introducti on 

Developments  in  high-energy  lasers  sometimes  come  faster  than 
the  reaction-rate-constant  experimental  data  necessary  to  analyre 
these  lasers.  For  example,  both  the  KrF  and  HgCl  lasers  are  candidates 
for  the  Department  of  Defense  High-Energy  Laser  Program.  Yet , no 
complete  analyses  of  their  reaction  kinetics  exist  to  determine  their 
scalability  to  larger  than  a laboratory  model.  These  analyses  do 
not  exist,  because  experimental  data  is  not  complete  on  all  the 
major  reactions  to  determine  their  reaction-rate  constants. 

As  more  potential  high-energy  lasers  are  developed,  the  need 
to  determine  reaction-rate  constants  increases.  Two  approaches  can  be 
taken  to  solve  this  expanding  problem.  One  method  to  determine  which 
react  ion-rat e constants  are  most  important  to  the  reaction  kinetics 
of  a laser  is  to  experimentally  determine  all  the  possible  reaction- 
rate  constants.  \ second  method  is  to  determine  theoretically  all 
the  contributing  reaction-rate  constants.  Neither  method  is  totally 
satisfactorv,  because  of  the  time  that  is  required  to  do  each  in 
detai 1 . 

A better  method  to  determine  the  important  reaction-rate  constants 
is  to  do  simple  theoretical  calculations  to  determine  an  order-of- 
magnitude  for  the  reaction  rates.  After  the  important  contributing 
reactions  to  the  laser  are  theoretically  determined,  experiments  can 

1 


be  conducted  to  ventv  that  these  reactions  are  important  and  to 


determine  the  exact  reaction-rate  constants  for  these  reactions.  In 

this  paper,  simple  theories  and  theoretical  calculations  are  reported 

for  the  HgCl  laser  reaction  kinetics. 

Both  Parks  (Ref  11  and  Lden  (Ref  21  report  observing  lasing  in 

0 

a mercury-chloride  laser  at  a wavelength  of  55"6  A.  The  lasing 


is  achieved  by  electron-beam  pumping  oi  mixtures  of  Ar,  Xe,  Hg, 


and 


Eden  reports  that  double  line  laser  oscillation  is  achieved 


at  55"6.2  A and  3583.5  A with  a mixture  of  S5 


An  oscillogram  of  one  of  Eden's  HgCl  laser 


actions  with  a single  laser  oscillation  at  55“6.2  A is 


l ven 


and  is  shown  in  Figure  1 


Fig.  1.  Oscillogram,  of  typical  HgCl(B-X)  fluorescence  (upper  trace! 
and  laser  emission  waveforms  for  a mixture  of  $5.4'’  Ar/10.4°$  Xe 
d-O^  Hg/0.2°i  CCl^  , pressure  = 3.43  atmospheres,  T,ejj  = 304°C,  and 

T . = 260°C.  Horizontal  scale:  20  ns/div;  vertical  units 

reservoir 

arbitrary  (Ref  2)  . 


Three  characteristics  of  the  HgCl  laser  make  it  an  excellent 
candidate  for  a high-energy  laser.  First,  the  wavelength  of  the  HgCl 


laser  is  55“e>  A which  is  in  the  blue-green  portion  of  the  spectrum. 

Laser  light  in  the  visible  range  is  attributed  with  better  atmospheric 
propagation  characteris t ics  than  laser  action  with  a longer  wavelength 
such  as  the  CO,  laser  at  10.6  microns.  Problems  such  as  atmospheric 
absorption,  aerosol  absorption  and  scattering,  thermal  blooming,  and  wind 
turbulence  are  greater  limiting  factors  on  lasers  with  a longer  wave- 
length (Ref  5).  Because  it  lases  in  the  visible  part  of  the  spectrum,  the 
capability  exists  for  the  HgCl  laser  to  be  used  during  a greater  variety  of 
atmospheric  conditions.  Second,  the  potential  quantum  efficiency  of  the 
HgCl  laser  is  45%  if  it  is  pumped  by  a discharge  (Ref  11.  Obviously,  this 
high  quantum  efficiency  makes  the  entire  system  more  efficient.  Third, 
the  energy  per  photon  of  the  HgCl  laser  is  higher  than  for  a laser  with 
a longer  wavelength. 

The  HgCl  potential  energy  curves  are  shown  in  Figure  2. 

Before  the  HgCl  laser  can  be  exploited  as  a high-eneigy  laser,  a 
detailed  analysis  of  the  HgCl  laser  reaction  kinetics  must  bo  completed. 

The  theories  and  calculated  theoretical  values  that  are  reported  in  this 
paper  provide  some  information  toward  determining  which  reactions  are 
important  in  the  HgCl  laser. 

The  approach  taken  to  analyze  the  HgCl  laser  in  this  paper 
is  a three  step  approach.  First,  a simple  model  of  the  reaction 
kinetics  for  the  HgCl  laser  is  constructed  using  Fden's  analysis 
(Ref  2)  as  a guide.  Second,  simple  theories  are  investigated 
for  the  various  reactions,  and  theoretical  calculations  are  made 
to  determine  the  important  reaction-rate  constants.  These  theoretical 


calculations  are  compared  to  available  experimental  data  to  determine 
the  validity  of  the  theories.  Third,  the  rate  equations  of  the  model 
are  solved  and  subjected  to  a sensitivity  analysis  to  further  determine 
which  reaction-rate  constants  need  more  research.  During  the  sensi- 
tivity analysis,  reaction-rate  constants  are  varied  to  determine  the 
sensitivity  of  the  entire  model  to  a given  reaction-rate  constant. 

The  investigation  of  reaction-rate  constants  is  limited  to  those 
applicable  to  the  llgCl  laser  model.  The  following  type  of  reactions 
are  included: 

1.  ionisation  by  electron  impact 

2.  electron-ion  dissociative  recombination 

3.  electron-molecule  dissociative  attachment 

4.  ion-neutral  association  reactions 

5.  bimolecular  charge  transfer  reactions 

o.  ion- ion  recombination 

7.  "harpoon"  reactions 

3.  excitation  by  electrons 

9.  radiative  decay  of  excited  atoms 

10.  Cl-Cl  recombination 

11.  quenching  of  excited  ligCl  by  Cl; 

The  following  general  assumptions  are  made.  First,  classical 
mechanics  and  quantum  mechanics  are  valid.  Second,  colliding  heavy 
particles  follow  the  trajectories  predicted  by  classical  mechanics. 

Third,  the  energy  transitions  of  the  internal  states  of  heavy 
particles  calculated  using  quantum  mechanical  perturbation  theory 
are  valid.  Heavy  particles  are  particles  other  than  electrons.  Fourth, 
the  perturbed  Hamiltonian  determined  by  the  time-dependent  collisional 
interaction  using  a classical  collision  model  applies  to  this  problem. 
Fifth,  the  electron  energy  distribution  and  the  gas  energy  distribution 
are  Maxwellian  and  are  at  different  average  energies  of  4eV  (Ret'  31 
and  0.0639eV,  respectively.  The  assumption  that  the  electrons  are  in 
a Maxwellian  distribution  is  made  to  simplify  calculations.  The  validity 


1 1 . Mode  1 


Intrcdu  ct  ion 

In  this  chapter,  a simple  model  is  developed  for  the  HgCl 
laser.  This  simple  model  is  used  as  a guide  for  investigation  of 
specific  reactions  in  the  HgCl  laser.  The  first  step  in  the  development 
of  this  model  is  a brief  discussion  of  Eden's  analysis  ( Re f 2).  The 
second  step  is  the  selection  of  the  reactions  that  are  included  in 
the  model.  Scire  of  the  reactions  that  are  not  included  in  this 
model,  yet  can  affect  the  lasing  action,  are  listed  in  Appendix  A. 

The  third  step  is  to  list  the  rate  or  continuity  equations  for  each 
species  that  is  involved  in  this  model.  Some  concluding  remarks  about 
the  model  are  included  as  the  last  section  of  this  chapter. 

Eden's  Analysis  (Ref  2) 

Eden  reports  that  his  experiments  were  conducted  primarily  at 
pressures  of  2 to  5 atmospheres.  At  these  high  pressures,  Eden  states 
that  harpoon  reaction  in  equation  1 is  not  significant.  The  harpoon 
reaction  creating  the  HgCl*  in  reaction  1 is  the  result  of  outer  shell 
electrons  rearranging  to  form  an  attractive  bond  between  the  Hg*  atom 
and  the  Cl  atom. 


Hg*(3P2)  + CC1U  -*  HgCl*  + CC13 


(1) 


+ Hg 


(2) 


Xe*(3P) 


Hg* ( 3S1  + Xe 


The  small  cross  section  for  reaction  2,  a * 10  5 A2  limits  the  Hg* 
that  is  created  to  insignificant  number  densities.  This  lack  of 
Hg*  necessary  for  reaction  I makes  this  mode  of  HgCl*  formation  insigni- 
ficant. Therefore,  the  primary  mode  for  HgCl*  formation  is  the  ionic 
recombination  of  Hg+  and  Cl*  by  reaction  5. 


Hg  * Cl  * M - HgCl* (B2E)  * M(M  = Ar.Xe) 


According  to  Eden,  the  Hg  is  produced  in  the  following  manner. 

The  electron  beam  energy  is  absorbed  by  the  Ar  to  create  Ar+  by  the 
electron-impact  ionization  of  Ar.  This  reaction  is  shown  in  reaction  4. 

P 

A high  energy  beam  electron  or  primary  electron  is  indicated  by  e 
An  electron  other  than  a beam  electron  is  indicated  by  e and  is 
referred  to  as  a secondarv  electron. 


p + P 

e + Ar  -*■  Ar  + e + e 


The  Ar  reacts  with  two  other  Ar  atoms  in  an  ion-neutral  recombination 
reaction  to  get  Art  , and  is  shown  in  reaction  5. 


Ar  + 2Ar  -*■  Ar?  + Ar 


The  Arj  then  transfers  a positive  charge  to  a Xe  atom  in  the 
bimolecular  charge  transfer  reaction  that  is  shown  in  reaction  6, 


Ar2  * Xe  Xe  ♦ 2Ar 


3 


The  excitation  of  the  metastable  Ar 


by  the  two-body  reaction  11. 


* is  transferred  to  a Xe  atom 


Ar*  + Xe  -*  Xe*  + Ar 


(11) 


The  metastable  Xe*  then  creates  excited  XeCl*  by  the  harpoon- like 
reaction  12. 


Xe*  * CClu  - XeCl*  * CC13 


(12) 


An  alternate  method  of  formation  for  XeCl*  is  the  ionic  recombination 
ot'  and  Cl  that  is  shown  in  reaction  13. 

* Cl'  + M - XeCl*  * M (M  = Ar.Xe)  (13) 


This  reaction  scheme  of  contributing  and  competing  reactions 
indicates  the  need  to  know  the  reaction-rate  constants  for  all  the 
reactions  does  exist.  A simple  analytical  analysis  of  all  possible 
reactions  is  impossible  in  a time-dependent  model.  Therefore,  a 
simple  model  is  constructed  to  analyte  some  of  the  reactions  that 
are  given  in  Eden’s  analysis  (Ref  2).  Some  reactions  that  are  not 
considered  by  Eden  are  included  in  the  model  for  analysis. 

Simple  Model 

The  model  constructed  in  this  section  is  developed  by  using 
Eden's  analysis  (Ref  2)  as  a starting  point.  Several  reactions  are 
added  to  Eden's  analysis  (Ref  2)  to  emphasite  specific  points  concerning 

10 


the  HgCl  laser.  Most  of  Eden's  analysis  for  high  pressures  (Ref  2) 
is  included.  Eden's  source  of  Xe*  is  modified  and  explained  later. 
Therefore,  reactions  5 through  13  are  included  in  this  model  with  the 
exception  of  reaction  11. 

In  this  model,  reactions  9 and  12  are  modified  from  Eden's 
analysis.  Eden  (Ref  2)  and  Parks  (Ref  1)  state  that  their  source 
of  Cl  was  the  electron-molecule  dissociative  attachment  of  CC1U 
shown  in  reaction  9.  This  model  uses  the  same  type  of  electron-molecule 
dissociative  attachment  reaction  to  create  Cl"  . The  source  molecule 
chosen  for  this  model  is  changed  from  CC14  that  is  used  by  Parks 

O 

and  Eden  to  Cl2  , because  the  cross  section  of  Cl^  is  90  A2 
(Ref  6)  for  the  harpoon  reaction  that  is  investigated.  The  cross 

O 

section  for  the  CC14  reaction  is  34  A2  (Ref  6).  The  Cl2  molecule 
is  also  substituted  for  CC1U  in  reaction  12.  This  slight  change  does 
not  present  any  theoretical  problems,  but  this  change  may  create 
experimental  problems  due  to  the  health  hazards  and  corrosive  nature 
of  Cl,  . Eden's  comment  concerning  the  slow  attacking  on  his 

cavity  windows  (Ref  21  indicates  a potential  problem  with  3ny  Cl 
atom  creating  chain.  The  modified  reactions  9 and  12  are  numbered 
9a  and  12a,  respectively. 

e ♦ Cl2  -*■  Cl"  + Cl  (9a) 

Xe*  * Cl2  - XeCl*  + Cl  (12a) 

In  this  model,  HgCl*  is  created  by  the  ion-ion  recombination 
of  Mg  and  Cl  in  reaction  3 and  the  harpooning  reaction  with 
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Hg*  and  Cl2  in  reaction  14  (Ref  6)  . 


Hg*(3P,)  + Cl2  - HgCl*  (B)  + Cl  (14) 

The  Hg*  necessary  for  reaction  14  is  not  created  by  reaction  2 
because  of  the  small  cross  section  for  transfer  of  the  Xe*  excitation 
to  Hg*  in  reaction  2 (Ref  2) . The  alternate  method  that  is  used  is 
the  electronic  excitation  of  Hg  in  reaction  15. 

e + Hg  -*  Hg*  (3P2)  + e (15) 

Only  Hg*  (3P2)  is  considered,  because  Hg*  (3P-)  is  found 
experimentally  to  be  the  dominant  mercury  atom  to  react  with  Cl2 
The  Hg*  ( - Pg ) reaction  cross  section  for  a harpoon  reaction  with 
Cl ^ to  the  IlgCl*(B)  state  is  found  to  be  one-fifth  the  Hg*  (3P2) 
reactive  cross  section  for  the  same  reaction.  Hg*  (3Pj)  does  not 
appear  to  create  a chemiluminescence  with  Cl2  (Ref  6). 

The  possible  loss  of  Hg*  (3P2)  by  electron-impact  ionization 
is  also  included  in  this  model.  This  reaction  is  a source  of  ilg+ 
and  electrons  as  is  seen  in  reaction  16. 

e + Hg*  (3P2)  - Hg+  - 2e  (16) 

For  the  remaining  text,  Hg*  is  used  to  symbolize  Hg*  (3P2)  unless 
a distinction  must  be  made  between  the  various  excited  states  of  mercury. 

12 


. - ; 


Sources  of  A r*  , Xe*  , Ar  , and  Xe  similar  to  those 
* ♦ 

listed  for  Hg  and  Hg  are  considered  in  this  model.  The  following 
reactions  of  Ar  and  Xe  with  secondary  electrons  are  considered. 


e «■  Ar  -*■  Ar*  + e (17) 

e + Ar*  -*■  .\r*  * 2e  (18) 

e + Xe  -*■  Xe*  + e (19) 

e ♦ Xe*  -*•  Xe+  + 2e  (20) 


A source  of  Xe*  not  mentioned  by  Eden  is  the  electron-molecule 
dissociative  recombination  of  Xe2  • This  reaction  is  theoretically 
the  same  as  the  electron-molecule  dissociative  recombination  of  Ar2 
in  reaction  10.  Reaction  21  shows  the  electron-molecule  dissociative 
recombination  of  Xet 

e + Xe£  -*•  Xe*  * Xe  (21) 

The  Ar*  and  Xe*  provide  competition  with  Hg*  for  use  of 
CI2  • The  formation  of  HgCl*  and  XeCl*  is  shown  in  reactions  14 
and  12a,  respectively.  The  reaction  of  Ar*  with  CI2  to  form 
ArCl*  (Ref  7)  is  shown  in  reaction  22. 

Ar*  ♦ Cl2  * ArCl*  ♦ Cl  (22) 
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Possible  effects  of  the  radiative  decay  of  XrCl 


XeCl 


and  HgCl*  are  included  in  this  model.  These  reactions  are  shown  in 
reactions  23,  24,  and  25  (Refs  2,  ",  3). 


ArCl*  - Ar  ♦ Cl  + hv  (1750  A) 


XeCl*  - Xe  + Cl  + hv  (5080  A) 


HgCl*  -*  HgCl  + hv 


Reactions  25  and  24  are  bound-free  type  transitions  (Ref  7) . 

The  existence  of  Ar*  and  Xe*  created  in  several  reactions  indicates 
possible  competition  with  Hg  for  neutralizing  reactions  with  Cl" 
atoms  created  in  reaction  9a.  Reaction  13  shows  the  Xe*  and  Cl" 
ion-ion  recombination  in  reaction  5.  Reaction  26  shows  Ar*  and  Cl" 
ion-ion  recombination  and  competes  with  reaction  3. 


Ar  + Cl  + M -*•  ArC  1 * + M (M  = Ar.Xe) 


Several  reactions  create  a neutral  Cl  atom.  Two  reactions  for  the 
three-body  neutral-neutral  recombination  of  Cl  to  Cl^  are  included 
in  the  model  to  investigate  the  impact  they  have  on  the  overall  laser 
action.  Reactions  2~  and  23,  which  were  investigated  by  Clyne  and 
Stedman  (Ref  9)  are  shown  here. 


atm m 


Cl  + Cl  + Ar  •*•  Cl-  + Ar 


Cl  + Cl  + Cl2  - 2C 1 2 


Eden  (Ref  2)  shows  the  primary  electrons,  e , for  the 
electron  beam  only  reacting  with  Ar  in  the  creation  of  Ar+ 
in  reaction  4.  Three  other  reactions  using  the  primary  electrons 
are  included  in  this  model.  The  excitation  of  Ar  by  a primary 
electron  is  shown  in  reaction  29.  The  excitation  of  Xe  by 
a primary  electron  is  shown  in  reaction  30.  Finally,  the 
ionication  of  Xe  by  a p-imary  electron  is  shown  in  reaction  51. 


I 


S 


p 

e + 


Ar  - Ar* 


P 

e 


(29) 


P 

e 


Xe  - Xe* 


P 

c 


(30) 


P 

e * 


Xe 


Xe^ 


P 

e 


♦ 


e 


(31) 


The  symbol  e represents  the  ionized  electron  and  is  referred  to  as 
a secondary  electron. 

The  final  reaction  that  is  included  in  this  model  is  a quenching 
reaction  of  HgCl*  by  Cl2  . This  reaction  is  shown  in  reaction  32. 

In  this  reaction,  the  HgCl*  is  quenched  to  a bound 

HgCl*  + Cl2  - HgCl  * Cl2  (32) 


state  for  HgCl  (Ref  S) . 

Several  facts  are  interesting  about  this  model.  First,  the  model 
does  not  include  the  radiative  decay  of  Hg*(3P2)  . This  radiative 

lifetime  is  extremely  long  and  results  in  a transition  probability 
of  10  1 sec  1 (Ref  10)  which  should  not  affect  this  problem.  Second, 
the  model  does  not  include  the  possible  existence  of  a stabilized 
molecule  of  Cl2  at  higher  operating  pressures  (Ref  11).  This 
molecular  ion  is  a possible  source  of  several  types  of  reactions. 

Third,  the  model  does  not  include  a loss  term  for  HgCl  molecules. 

This  fact  creates  some  miner  problems  in  the  analytical  equilibrium 
soluti-en  of  the  model.  Essentially,  no  true  equilibrium  solution  exists 
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for  this  simple  model,  since  at  least  one  species  is  always  increasing 
or  always  decreasing  in  number  density.  Appendix  3 explains  the  lack 
of  photodetachment  of  Cl" 

As  a summary,  the  following  reactions  are  included  in  this  simple 
model  of  the  HgCl  laser. 


P + p 

e ♦ Ar  -*  Ar  + e + e (4) 


P P 

e + Ar  -*■  Ar*  + e 


P P 

e + Xe  Xe*  + e 


P V + P 

e + Xe  ■*  Xe  + e + e 


Ar  + 2Ar  -*■  Ar,  + Ar 


Ar?  + Xe  - Xe  + 2Ar 


Xe  + Xe  + Ar  -*•  Xe-  + Ar 


Xe  + Xe  Xe  - Xe?  + Xe 


Xej  + Hg  -*•  Hg  + 2Xe 


e + C1t  -*  Cl  -►  Cl 


(291 


(30) 


(31) 


(5) 


(6) 


(?a) 


(“bl 


(S) 


(93) 


r 
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Ar*  + Cl'  + M -*■  ArC  1 * «■  M (26) 

ArCl*  ■*  Ar  Cl  hv  (23) 

XeCl*  -*•  Xe  * Cl  * hv  (24) 

HgCl*  -v  HgCl  ♦ hv  (25) 

Cl  t Cl  ♦ Cl2  * 2C12  (2S) 

Cl  + Cl  + Ar  -*•  Cl2  * Ar  (27) 

HgCl*  + Cl2  -*  HgCl  + Cl2  (32) 


Rate  Equations 

Each  reaction  in  the  model  occurs  with  a specific  reaction-rate 
constant.  This  reaction-rate  constant  multiplied  by  the  number 
densities  of  the  reacting  species  equals  the  number  per  cubic  centimeter 
per  second  of  a given  species  that  is  lost  or  created  in  a given  reaction. 
Rate  equations  must  be  determined  for  each  species  in  the  model,  so  the 

interdependence  of  all  the  species  can  be  resolved  in  the  final  solution. 

» ;* 

In  this  section,  rate  equations  are  listed  for  all  species  included 
in  the  model.  A specific  species  is  any  atom,  molecule,  or  electron 
with  a specific  chemical  characteristic.  For  example,  Hg*  , Hg  , 
and  IIg*(3P~)  are  all  different  species.  The  reaction-rate  constants 
in  this  paper  are  symbolized  by  , where  i is  the  reaction  number 
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used  in  the  previous  sections.  For  example,  the  reaction-rate  constant 
for  reaction  4 is  R4  . The  number  density  for  each  species  is  indicated 


dv  n , 


where  A is  a specific  species.  For  example,  the  number 


densitv  of  Ik*  is  indicated  bv  n,,  * . The  time-rate  of  change  of 

Hg 


a specific  species  is  indicated  by  n 


A 


where  n indicates  the 


different ial 


dn 


and  A indicates  a specific  species.  For  example, 

the  time-rate  of  change  of  Cl"  is  indicated  bv  n , 

* -cl* 


Other  svmbols  introduced  in  this  section  are  S 


eh 


eb 


% 


.Ar 


\e 


Ar 


and 


Vr+ 


*Xe 


eb 

The 


symbol  + is  the  number  of  Ar  ions  or  secondary  electrons  created 
per  cubic  centimeter  per  second  by  the  incident  electron  beam  in 
reaction  4. 


- . •’eb  Mr 
Ar  " q ‘\\r 


(531 


The  svmbol  i . indicates  the  current  of  the  incident  electron  beam. 
J eb 

The  symbol  n^  represents  the  number  density  of  beam  electrons  in 
the  electron  beam.  The  symbol  v ^ represents  the  speed  of  a 


monoenergetic  primary  electron. 


j , = q n . v , 

eb  le  eb  eb 


(34) 


Substitution  of  equation  34  into  equation  35  results  in  equation  33a. 


S,  + = n,n,  v , o, 

Ar  eo  Ar  eb  Ar 


(33a  1 
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Equation  13a  is  in  a more  standard  rate  of  reaction  form.  The  number 

densities  of  reacting  species  are  n ^ and  n^r  . The  reaction-rate 

constant  is  j,  v . where  j,  is  the  cross  section  for  the 

ionitation  of  Ar  by  primary  electrons  with  a specific  energy  and 

is  the  speed  of  the  approaching  electron.  In  this  case,  v ^ is 

approximately  the  speed  of  approach  of  the  two  particles,  because  the 

speed  of  the  electron  is  several  orders -of -magnitude  greater  than  the 

speed  of  the  Ar  . Similar  arguments  exist  for  the  terms  including 

jv  , , and  <j..  where  jv  is  the  cross  section  for 

ionication  of  Xe  by  a primary  electron  in  reaction  51,  is  the 

\e 

cross  section  for  excitation  of  Xe  to  a metastable  state  by  a primary 
electron  in  reaction  50,  and  j r is  ti'e  cross  section  for  excitation 
of  Ar  to  a metastable  state  by  a primary  electron  in  reaction  29. 

The  rate  equations  for  secondary  electrons,  Ar,  and  Xc  are 
presented  in  equations  55,  56,  and  57,  respectively. 


n = v , a.  n , n.  v , <jv  n , nv  - R,.  n n . 

e eb  Ar  eb  Ar  eb  Xe  eb  Xe  9a  e cl- 


+ R,,  n n..  * + Rl0  n n. 

lb  e Hg*  13  e Ar’ 


R,n  n nv  * 

20  e Xe*  (oa) 


R,.  n n,  + 
10  e Ar, 


R„  n nv  + 
21  e \e-> 


n.  = -v  . , n , n,  - v . o. , n , n.  - R_  n + nr 

Ar  eb  Ar  eb  Ar  eb  *Ar  eb  Ar  s Ar  Ar 


(.06 1 


+ 2R,  n.  + n,,  - R,_  n n.„  + R,„  n n + + R,_  n. 

b Ar?  Xe  1/  e Ar  10  e Ar^  2o  ArCl* 
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v , ov  n n.. 
eb  Xe  eb  Xe 


Xe 


v , o*v  n , nv 
eb  *Xe  eo  Xe 


R,  n,  + nv 
6 Ar->  Xe 


(R_  n,  + R_,  n„  ) nv,  n..  * * 2R„  nv  + n 

a Ar  b Xe  Xe  Xe  8 Xe-. 


Hg 


(3~) 


- R 
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n 

e 


n 


Xe 


n 

e 


n 


XeC  1 


■k 


The  term  (R_  n,  + R„u  nv  ) nv  n.,  + that  is  found  in  the  rate 
equations  for  Xe  , Xe+  , and  Xet  is  the  sum  of  the  two  possible 
contributions  from  reactions  "a  and  7b  to  these  specific  rate  equations. 
This  term  is  the  result  of  two  possible  third  reacting  bodies  in  the 

formation  of  Xe-  . The  third  body  is  not  included  in  the  ionic 

. + +■  + » , 
recombination  ot  Ar  , Xe  , and  Hg  with  Cl  tor  reasons 

that  are  discussed  in  Chapter  3. 

The  rate  equations  for  the  remaining  species  involved  in  the 

model  are  listed  in  Equations  58  through  53. 


n.,  + = v . n . nv  + R.  n .♦  n,.  - (R,  n * R_,n..  )nv  nv  * 

Xe  eb  Xe  eb  X e 6 Ar?  Xe  , a \r  b Xe  Xe  Xe 


* R,_.  n nv  * - R, . nv  + n • 

20  e Xe*  lo  Xe  cl 


(SSI 


n y + 
Xe-. 


(R„  n,  * R_,  nv  In..  nv  « 
7a  Ar  b Xe  Xe  Xe 


R_  n . + n..  - R,,  n nv  * 

S Xe-  Hg  21  e Xe- 


(39) 


n,  + = v . a. 

Ar  eb  Ar 


n . n.  - R_  n,  ♦ n.z  * R,.,  n n - R n * n - f 101 

eb  Ar  5 Ar  Ar  15  e Ar*  25  Ar  cl 


n + = R.n+n,- 


Ar 


? 3 


Ar  Ar 


R . n,  n„  _ R,  , n n,  + 
6 Ar-,  Xe  ' 10  e Ar-. 


(41) 


n,  * = v ,q..  n ,ii,  + R,„n  n,  - R,0n  n,  , * Rlrtn  n,  * 

Ar*  eb  *Ar  cb  Ar  17  e Ar  18  e Ar*  10  e Ar-. 


R..n,  „n„, 
22  Ar*  Cl 


142) 


n,,  „ = v uO*v  n ,n,.  + R,ni\  nv  - R,,,n  n„  . + R,,n  n + 

\e*  eb  *Xe  eb  \e  19  e Xe  20  e Xe*  21  e Xe? 


R12a\e*nCl, 


(43) 


Hg 


-R0nv  +n,, 

8 Xe2  Hg 


R,.n  n,. 
Is  e Hg 


44) 


nHg+ 


R.n,.  +nn  + R-w-n  n,,  . 
S \e?  Hg  16  e Hg* 


R_n  +n  • 
->  Hg  c 1 


(45) 


Hg' 


R,  _n  nLI 
la  e Hg 


R,,n  n..  * 
16  e Hg 


R , • n,  * ii/-« , 
14  Hg*  Cl; 


(46) 


Cl, 


'R9anenCl?  + R28nClnCl,  * R27nClnAr 


(47) 


n„,-  = R„  n n 


Cl  9a  e Cl2  3 Hg 


- R.ri,,  .n„,_  - R,  ,nv  +n_,-  - R,,n,  +n„ 

, a Hg+  Cl  la  Xe  c.l  2b  Ar  u 


Cl 


(48) 


"«  * V."ci,  * !uV"ci:  * R::nAr-"a2  * #i:a"xe-”ci2 


' XV'vCl*  * R24nXeCI*  ' R:SnClnCl2  ' !*;7nCl"Ar 


n,.  r1,  = R,  inu,*n^,  + R,n..  +n,  - - R.^n.,  r,+nri  - R,rnu  ( 

HgC 1*  14  Hg*  Cl2  o Hg  LI  o2  HgCl*  Cl?  25  HgCl* 


"\rCl*  = R22nAr*nCl,  + R26nAr+nC 


Cl 


R25nArCl* 


nv  = R,,  nv  *n_,  + R,_nv  +iu,-  - R,,nv. 

Xetl*  12a  Xe*  Ll?  lo  Xe  Cl  24  Xetl* 


nHgCl  = R32nHgCl*nCl2  + R25nHgCr 


Conclusion 

The  simple  model  that  is  presented  here  does  not  include  all  the 
possible  reactions  taking  place  in  the  HgCl  laser.  This  model  does 
include  several  reactions  not  considered  by  Parks  (Ref  1)  and  Eden 
(Ref  2).  A major  deviation  from  Parks'  and  Eden's  experiments  is  the 
substitution  of  Cl?  for  CC14  in  the  model. 

Before  any  solution  can  be  found  for  the  rate  equations  that  are 
presented  in  this  section,  the  reaction-rate  constants  for  each  reaction 
must  be  found.  The  reaction-rate  constants  are  developed  in  Chapter  3, 
and  a solution  for  the  rate  equations  is  presented  in  Chapter  4. 
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III.  Reaction-Rate  Constants  and  Theories 

Overview 

In  this  chapter,  the  reaction-rate  constants  necessary  to  solve 

the  rate  equations  in  Chapter  2 are  developed.  The  processes  that 

are  discussed  and  include  theoretical  sections  are  dissociative 

recombination,  dissociative  attachment,  ion-neutral  association 

reactions,  charge  transfer  reactions,  ion- ion  recombination,  and 

harpooning  reactions.  Reaction-rate  constants  that  are  presented 

based  on  experimental  data  are  secondary-electron  excitation  of 

Ar  , Xe  , and  Hg  , ionization  of  excited  states  by  secondary- 

electron  impact.  Cl -Cl  recombination  reactions,  radiative  decay  of 

ArCl*  , XeCl*  , and  HgCl*  , ionization  and  excitation  of  Ar 

and  Xe  by  the  electron  beam,  and  HgCl*  quenching  by  Cl- 

Some  of  the  reaction-rate  constants  for  these  processes  are 

the  result  of  direct  measurement  of  the  reaction-rate  constant,  and 

others  are  the  result  of  cross  section  measurements . If  only  cross 

sections  are  found,  the  reaction-rate  constant  for  that  process  can 

be  calculated  by  multiplying  the  cross  section  times  the  average 

relative  velocity  of  the  two  particles.  The  average  relative  velocity, 

v , is  found  for  oarticles  in  a Maxwell -Boltzmann  distribution  to 
r 

3KT  * /o 

be  2 where  k is  the  Boltzmann  constant,  T is  the  temperature 

r 

of  the  gas,  and  M is  the  reduced  mass  of  the  particles. 


v- 


Li A 


1",  13)  exist  for  dissociative  recombination  reactions  and  orovide 


reasonable  data  for  reaction-rate  constants.  This  section  includes 
summaries  of  applicable  theoretical  and  experimental  work  that 
exist  for  dissociative  recombination. 

Theory.  Dissociative  recombinat ion  requires  a free  electron 
with  positive  energy  to  give  away  energy  to  become  a bound  electron 
with  negative  energy  (Ref  12).  This  excess  energy  can  be  radiated 
away  or  given  up  in  a collision  with  a third  body.  Dissociative 
recombination  is  essentially  a two-body  collision  between  an  electron 
and  a molecule  containing  two  atoms.  One  of  these  atoms  is  the 
third  body  necessary  to  absorb  the  excess  energy  of  the  electron. 

The  electrons  that  are  captured  in  dissociative  recombination 
are  thermal  electrons  (Ref  12).  These  electrons  provide  the  energy 
necessary  for  molecular  ions  to  form  unstable  intermediate  neutral 
molecules  AB"  that  dissociate  into  two  neutral  atoms. 

Figure  3 shows  the  shape  of  potential  energy  curves  necessary 
for  direct  dissociative  recombination  to  occur.  The  internuclear 
separation  distance  is  the  point  at  which  the  energy  the  thermal 

electron  provides  for  the  reaction  is  sufficient  to  move  the  bound 
molecular  ion  to  the  higher  potential  energy  curve  of  the  unstable 
neutral  molecule  AB"  . This  point,  , is  called  the 

capture  point. 


1 r 
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In  this  chapter,  the  electrons  are  still  assumed  to  be  in  a 
Maxwell-Boitzmann  distribution  with  an  average  energy  of  4 kTg 
equal  to  4 eV  . This  assumption  is  made  to  simplify  calculations. 

A solution  of  the  Boltzmann  equation  is  required  to  find  the  actual 
distribution  function  of  electrons.  This  new  distribution  is  required 
to  give  correct  answers  for  reaction-rate  constants  that  are  calculated 
for  electron- impact  excitation  and  ionization  by  secondary  electrons. 

A solution  of  the  Boltzmann  equation  is  beyond  the  scope  of  this  paper. 

The  next  section  of  this  paper  contains  a discussion  of 
dissociative  recombination.  The  reaction-rate  constants  for 
dissociative  recombination  are  experimentally  determined. 

Dissoc iative  Recombination 

Introduction . Dissociative  recombination  is  the  reaction  n 
which  a diatomic  molecular  ion  captures  an  electron  and  dissociates 
into  two  neutral  atoms.  This  reaction  is  shown  in  reaction  54  where 
A and  B are  atoms. 


e ♦ AB  •*  AB  ' A + B (541 

The  AB"'  molecule  is  an  intermediate  unstable  molecule.  Two 
reactions  in  this  model  involve  dissociative  recombination.  In  these 
reactions,  Art  and  Xet  dissociate  to  form  the  appropriate  products 
shown  in  reaction  54. 

The  theoretical  work  (Refs  15,  13,  14,  15)  concerning  dissociative 
recombination  does  not  provide  simple  formulas  for  quick  calculation 
of  reaction-rate  constants.  Fortunately,  several  experiments  (Refs  16, 
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Fig.  3.  Schematic  representation  of  the  potential  energy  curves 
involved  in  the  direct  dissociative  recombination  process.  In  this 
simplified  example,  only  one  potential  curve  for  the  molecular  ion 
AB+  (with  an  electron  at  rest  at  infinity)  and  one  for  the  unstable 
molecule  AB"  are  shown.  Inset:  magnified  view  of  the  potential 
curves  in  the  curve-crossing  region.  (Ref  12). 


After  the  electron  is  attached  bv  the  molecular  ion,  two 


is  for  the  unstable  molecule  to  autoionice  and  eject  an  electron 


s the  unstable  molecule  back  into  the  bound  molecula 


into  two  neutral  atoms.  This  action  is  the  c 


recombination  reaction.  Figure  3 


of  the  molecular  ion  and  the  unstable  molecule  crossi 


point  R . This  ooint  is  called  the  stabil ization  noir.t  and  is 
the  internuclear  separation  distance  that  limits  the  time  tor  auto- 
ionization. After  the  neutral  atoms  are  separated  by  a distance 
greater  than  R_  , additional  energy  is  necessary  for  the  system 
of  particles  to  return  to  the  bound  molecular  ion  state  where  auto- 
ionization can  take  place.  Using  the  Born-Oppenheimer  approximation 
autoionization  is  impossible  at  internuclear  separations  greater 
than  . (.Ref  IT)  The  time  necessary  for  autoionization  is 

typically  10  14  sec,  and  the  time  necessary  for  the  neutral  atoms 

to  separate  to  greater  internuclear  separation  than  R.  is  10'  3 

1 5 

to  10  sec  (Ref  12).  Therefore,  the  dissociative  recombination 

process  is  more  likely  to  result  in  two  neutral  atoms  than  for 

the  unstable  molecule  to  autoionize. 

The  cross  section  for  the  dissociative  recombination  of  a 

molecular  ion  into  two  neutral  atoms  is  given  in  equation  55  where 

0(c)  is  the  cross  section  for  dissociative  recombination  as  a 

function  of  incident  electron  energv  t , o (c)  is  the  cross 

cap 

section  for  the  capture  of  an  electron  by  the  molecular  ion  to  form 
the  intermediate  unstable  molecule  AB"  , .and  S(s)  is  the 
probability  that  the  intermediate  unstable  molecule  survives  long 
enough  to  dissociate  and  is  called  the  survival  factor  vRef  IT). 


o(c)  = j (e)S(e) 

cap 


(55) 


Previous  arguments  indicate  that  5(c)  is  high  for  the7-mal  electrons 
and  can  be  considered  approximately  one. 


o 


The  3 tern  is  now  examined.  An  assumntion  is  made  bv 

cap 

Bards  ley  and  Biondi  that  the  unstable  molecule  AB"  is  formed  by 
the  excitation  of  a target  electron  to  create  a higher  excited  state 
and  the  incident  electron  falls  into  an  unoccupied  orbital  (Ref  12). 
issuming  Y(R)  represents  the  electronic  matrix  element  governing 


the  configuration  interaction,  ;(r,R)  represents  the  electronic 

wave  functions  of  the  molecules,  ;iR)  represents  the  nuclear  wave 

functions,  ?^(.r,R)  represents  the  incident  electron  wave  function, 

and  r represents  the  ratio  of  multiplicities  of  the  intermediate 

molecule  AB"  to  the  initial  molecular  ion  AB+  , 7 fc)  is 

cap^ 

given  in  equation  56. 


?capU)  = C2t3/W')(y)  !<:ab+(R)  !v(R)  ■ CAB  — (R): 


(56) 


The  ratio  of  multiplicities  is  divided  by  two  to  allow  for  the  two 
possible  spin  orientations  of  the  incident  electron  (Ref  12).  The 
electronic  matrix  element  Y(R)  is  given  in  equation  5“  where 
H ,(r,R)  is  the  Hamiltonian  for  the  molecular-electronic  wave  function, 

Cl 


'•'(R)  ■ <)AB+(r,R)JJr,RllHol(r,R)|^R„(r,Rl> 


el 


AB ' 


(57) 


The  (1/c)  proportional ity  for  7,  (e)  is  developed  in  reference  15. 

Figure  4 shows  examples  of  the  expected  nuclear  wave  functions 
for  AB+  and  AB"  . 
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Fig.  4.  Hypothetical  potential  energy  curves  and  associated  wave 
functions  for  one  of  the  states  of  the  molecular  ion  AB+  and  of  the 
unstable  molecule  AB"  involved  in  the  direct  dissociative  recom- 
bination process.  (Ref  12) 

The  ;vb_(R)  wave  function  shows  a large  probability  of  AB" 

particle  bing  near  the  classical  turning  point  R . This  wave 

C 

function  oscillates  rapidly  and  decreases  exponentially  is  the  inter- 

nuclear  separation  distance  R grows  greater  than  R_  . Bardslev 

h 

and  Biondi  (Ref  12)  approximate  ; ->(R)  with  a delta  function 

Ad 

located  at  R (Ref  19) . By  putting  the  delta  function  at  R , 

C ' c 

Bards  ley  and  Biondi  require  that  R^  be  close  to  R^  . The  delta 

function  is  given  in  equation  3$  where  U'  is  the  slope  at  R 


The  requirement  for  the  location  of  R^  is  given  in  equation  59 
where  U indicates  the  potential  energies  of  the  molecules. 


UAB~<Rc> 


AB' 


<RC> 


(59) 


The  ground-state  nuclear  wave  function,  ; + (R)  is  approximated 

Ad 

by  a harmoni  ilator  wave  function.  Substitution  of  these 

approximations  into  equation  56  yields  equation  60. 
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cap 


(e) 


(:*r3/f me)  (J)  (jpO  |V(R.)  |2  |c  VR+  (Rj|2 


(60) 


Bv  defining  the  capture  width  r as  (^-)  : V (R^J | 2 , equation  60 

reduces  to  equation  61. 
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This  development  indicates  that  R must  be  within  the  range  of 

c 

nuclear . vibrations  so  the  probability  of  capture  is  greatest. 

The  reaction-rate  constant  for  dissociative  recombination  is 
obtained  by  integrating  over  all  energies  for  the  product  of  the 
electron  distribution,  f(e)  , the  dissociative  recombination 
cross  section,  a (e)  , and  the  electron  velocity  as  'n  equation  62. 
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(62) 


If  f(e)  is  a Maxwel 1 -Boltzmann  distribution  of  electron  energies, 
equation  62  reduces  to  equation  63  (Ref  12). 
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Bardsley  and  Biondi  (Refs  12,  14,  15)  present  four  assumptions 
that  allow  the  dissociative  recombination  rate  for  the  ground 
vibrational  and  electronic  state  to  be  expressed  in  terms  of  electron 
temperature  . These  assumptions  are: 

(i)  The  survival  factor,  Sis)  , is  almost  one. 

(ii)  The  incident  electron  energy  is  less  than  1 eV. 

(iii)  U,0*-  is  such  that  R is  close  to  R_ 

AB  c E 

(iv)  The  vibrational  wave  function  for  AB+  can  be 

approximated  by  the  ground  state  harmonic  oscillator 

wave  function. 


The  resultant  dissociative-recombination  reaction-rate  constant 
is  given  in  equation  64  where  a is  the  amplitude  of  initial 
vibrational  motion,  R^  is  the  equilibrium  nuclear  separation 
distance,  and  T is  in  units  of  °K  (Ref  12). 


rF  (R  -R  )-  T 7 , 

W ■ s-su4e*p  <- aw- 

a" 


N'ot  all  experimental  data  agree  with  Te"5,5  dependence  in 
equation  64.  This  is  the  fact  that  led  Bardsley  (Ref  15)  to 
investigate  an  alternate  method  of  dissociative  recombination.  The 
alternate  method  is  called  indirect  dissociative  recombination  and 


is  represented  in  reaction  65. 


e + AB+  ->  AB'  -*■  AB"  - A + B 


(65) 


The  potential  energy  curves  representing  reaction  65  are  given  in 
Figure  5. 


Fig.  5.  Hypothetical  potential  energy  curves  for  the  molecular 
ion  AB+ , a highlying,  stable  Rydberg  state  AB'  , and  an  unstable 
molecule  curve  AB"  involved  in  the  indirect  dissociative  recombination 
process.  The  v = 0 vibrational  level  of  the  molecular  ion  AB+ 
and  the  v'  =,  0 and  1 vibrational  levels  of  the  Rydberg  state  AB  ' 
are  also  shown  (Ref  12) . 

In  Figure  5,  the  energy  contributed  bv  the  electron  is  =-  . This 

~n 

energy  must  be  close  to  the  energy  that  is  required  to  put  the  AB+ 
state  into  a AB'  Rydberg  state.  The  electron  energy  is  contributed 
to  the  nuclear  motion  of  the  AB'  state.  The  excited  Rydberg  state 
predissociates  to  AB"  because  of  the  location  of  the  potential 
energy  curve  crossing,  and  normal  dissociative  recombination  takes 


54 


place.  The  electron  temperature  dependence  for  the  reaction-rate 
constant  for  reaction  65  is  given  in  equation  66  where  n is  the 
possible  number  of  Rydberg  states. 


exp  (-e^kTJ 


(66) 


The  possibility  exists  for  direct  and  indirect  dissociative  recom- 
bination for  any  given  reaction.  Therefore,  the  eventual  theoretical 
temperature  dependence  can  vary  depending  on  the  available  Rydberg 
states  and  the  electron  temperature. 

Experimental  values  and  calculations.  O’Mallv  et  al.  (Ref  16) 
find  that,  if  the  location  of  the  stabilization  point  is  shifted 
slightly,  a better  probability  for  dissociative  recombination  exists 
if  the  ionic  molecule  is  in  an  excited  vibrational  state.  The  exper- 
imental data  for  neon  is  best  fit  by  a low  vibrational  state  model 
using  parameters  that  are  determined  by  a computer  program  (Ref  16). 

This  model  predicts  the  reaction-rate  constant  for  neon  to  vary 
according  to  T^  and  five  set  parameters. 

Fortunately,  experimental  data  for  the  dissociative  recombination 
of  Art  and  \e+  as  a function  of  incident  electron  energy  (Refs  17, 
18)  exists,  and  reasonable  predictions  for  these  reaction-rate  constants 
can  be  made  simply.  Mehr  and  Tiondi  (Ref  17)  find  that  the  reaction- 
rate  constant  for  Art  varies  as  Te-0-67  for  a gas  temperature  of 
500°  K,  Assuming  that  the  0.05  eV  increase  in  gas  temperature  of  this 
model  over  Mehr  and  Biondi's  (Ref  17)  experiment  does  not  affect  the 


vibrational  distribution  of  the  molecular  ions  greatly,  Meh-  and 


Biondi’s  T dependence  is  used  to  predict  a dissociative- 

recombination  reaction-rate  constant  of  3.9  x 10  ■ ctnVsec  for 

Art  in  this  model.  Shiu  et  al.  (Ref  IS)  find  that  for  Xe2  the 

dissociative-i'ecombination  reaction-rate  constant  varies  as  T ~3-’T- 

e 

for  T&  > 1000°  K and  a gas  temperature  of  300°  K.  Again,  assuming 

the  data  that  is  reported  for  a gas  temperature  of  300°  K are  valid 

for  the  temperature  of  the  model  of  333°  K,  a dissociative-recombination 

reaction-rate  constant  of  S.16  x 10  3 cm3/sec  is  calculated  for 
+ 

Xe2 

Conclusion . The  theory  for  dissociative  recombination  is  difficult 
to  use  for  predicting  reaction-rate  constants.  Fortunately,  sufficient 
experimental  data  exist  and  can  be  used  to  predict  dissociative- 
recombination  reaction-rate  constants  for  higher  electron  energies. 

The  data  used  for  these  predictions  are  for  electron  temperatures 
to  3000°  K for  Xe2  and  electron  temperatures  to  10000°  K for 
Art  (Refs  17,  13). 

The  assumption  that  data  for  gases  at  300°  X are  valid  for  gases 
at  533°  X is  probably  valid.  The  vibrational  spacings  of  atmospheric 
diatomic  ions  are  on  the  order  of  0.2  eV  (Ref  12).  The  energy 
difference  between  300°  K and  533°  K is  approximately  0.03  eV  and 
probably  is  not  near  a vibrational  level.  Therefore,  the  reaction-rate 
constants  for  the  dissociativerecombination  of  Art  and  Xet  are 
at  least  close  to  the  values  tnat  are  predicted  here. 

The  reaction-rate  constant  values  of  3.9  x 10  9 cm3 /sec  for 
Art  and  3.16  x 10  9 cm3 /sec  for  Xe,  are  used  in  Chapter  4 for 
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calculations  involving  an  applied  electric  field.  With  an  applied 
electric  field  across  the  plasma,  the  average  electron  temperature 
is  found  to  be  approximately  4 eV  (Ref  4) . For  electrons  at  these 
higher  energies  that  are  not  in  a Maxwell- Bo ltcmann  distribution,  it 
is  best  to  measure  the  reaction-rate  constant  rather  than  the  cross- 
section  for  the  reaction. 

For  calculations  in  Chapter  4 that  are  made  for  no  applied 
electric  field,  the  electrons  are  assumed  to  be  in  thermal  equilibrium 
with  the  gas.  The  dissociative-recombination  reaction-rate  constants 
that  are  used  in  this  case  are  5.99  x 10”  cm3/sec  for  Aro  and 
1.88  x 10  5 cnr/sec  for  Xe^  . This  assumption  is  discussed  in 
Chapter  4. 

The  next  section  is  about  a similar  process  to  dissociative 
recombination  that  is  called  dissociative  attachment. 

Dissociative  Attachment 

Introduction.  Dissociative  attachment,  a process  similar  to 
dissociative  recombination,  involves  the  capture  of  an  electron  by 
a neutral  molecule.  The  molecule  then  dissociates  to  form  a negative 
atom  and  a neutral  atom.  The  dissociative-attachment  reaction  in 
this  model  is  reaction  9a  shown  below. 


e + CI2  Cl  + Cl  (9a) 

The  theoretical  description  given  in  this  section  is  only  a 
general  presentation  that  is  based  on  Bards  ley  and  Mandl  (Ref  CO)  and 
the  potential  energy  curves  for  Clj  that  are  developed  by  Gilbert  and 

5? 
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Wahl  (Ref  21).  The  experimental  data  that  is  available  for  comparison 
to  theory  is  primarily  for  dissociative  attachment  of  Cl2  by  thermal 
electrons  (Refs  11,  22,  23,  24). 

Theoretical  and  experimental  results.  Dissociative  attachment 
is  the  capture  of  electrons  by  neutral  molecules  and  the  dissociation 
of  these  now  negative  molecular  ions  each  into  a negative  atom  and  a 
neutral  atom.  The  one  difference  between  dissociative  recombination 
and  dissociative  attachment  is  the  requirement  that  the  neutral  molecule 
involved  in  dissociative  attachment  must  contain  an  atom  with  a high 
electron  affinity.  Electron  affinity  is  the  energy  released  when  the 
neutral  atom  gains  an  electron.  Gilbert  and  Wahl  (Ref  21)  find  the 
electron  affinity  of  ground  state  Cl2  to  be  2.19  eV. 

Bardsley  and  Mandl  (Ref  20)  describe  dissociative  attachment  in 
a similar  manner  to  dissociative  recombination.  The  energy  of  the  cap- 
tured electron  must  be  great  enough  to  excite  the  target  molecule  to 
an  energy  greater  than  that  of  the  potential  energy  of  the  intermediate 
molecular  ion  as  is  shown  in  Figure  6 (Ref  20).  The  experimental 
results  of  Frost  and  McDowell  (Ref  25)  provide  some  justification  for 
this  theoretical  model.  The  experimental  results  give  a threshold 
electron  energy  for  the  dissociative  attachment  of  Cl2  to  be  1.6  eV 
and  support  the  supposition  that  this  dissociative  attachment  is  the 
result  of  formation  of  Cl2  (^ir  ) state. 

Recent  experiments  (Refs  11,  22,  23)  for  thermal  energy  electrons 
show  that  the  threshold  energy  for  the  dissociative  attachment  of 
Cl,  is  not  1.6  eV.  In  fact,  this  experimental  data  shows  that  thermal 
electrons  are  attached  to  CI2  rapidly.  Christodouiides  et  al.  (Ref  11) 
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Fig.  6.  Potential  energy  curves  and  nuclear  wave  functions  in 
dissociative  attachment.  The  final  wave  function  is  schematic 
only  (Ref  20) . 

find  the  dissociative-attachment  reaction-rate  constant  for  Cl,  to 

-in 

be  3.1  x 10  " cm-/sec.  Schultes  et  al.  (Ref  22)  report  a broad 

maximum  exists  for  the  reaction-rate  constant  for  dissociative 

attachment  of  Cl2  between  electron  energies  of  0.25  eV  and 

0.5  eV  . Sides  et  al.  (Ref  23)  report  a dissociative-attachment 

_a 

reaction-rate  constant  for  Cl,  of  (3.'T  z l.S)  x 10  ' cm3/sec  for  an 
electron  temperature  of  350°  K.  These  experiments  (Refs  11,  22,  23) 
confirm  the  formation  of  Cl,  (22+)  for  thermal  electrons. 
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The  potential  energy  curves  that  are  presented  by  Gilbert  and 
Wahl  (Ref  21)  are  shown  in  Figure  7.  These  curves  show  the 
intersection  of  the  Cl2  potential  energy  curve  by  the  C12(2I^) 
potential  energy  curve  at  the  minimum  of  the  Cl2  potential  well. 

The  possible  formation  of  Cl2  (2~*)  by  thermal  electron  excitation 
from  the  Cl2  potential  curve  to  the  C12(2Z*)  potential  energy 
curve  is  apparent.  The  possibility  of  excitation  from  the  Cl2 
potential  energy  curve  to  higher  levels  of  fl2  does  exist. 

Conclusion.  The  potential  energy  curves  for  Cl2  of  Gilbert 
and  Wahl  (Ref  21)  imply  that  a graph  of  the  dissociative-attachment 
cross  section  for  C1-,  versus  electron  energy  mav  contain  several 
peaks  corresponding  to  each  excited  potential  energy  curve  of  Cl2 
The  location  of  these  peaks  depends  on  the  electron  energy  required 
to  excite  the  Cl2  to  Cl2  . Further  research  is  required  to 
determine  the  electron  energy  dependence  of  the  dissociative-attachment 
reaction-rate  constant  for  Cl2 

The  reaction-rate  constant  for  dissociative-attachment  of  Cl 

2 

by  thermal  electrons  that  is  determined  by  Sides  et  al.  (Ref  23), 

3.7  x 10  ' cm3/sec  , is  used  in  this  model. 

The  next  section  deals  with  the  ion-neutral  association  reactions. 


Ion-Neutral  Association  Reactions 

Introduction.  In  this  section  two  theories  are  discussed  for 
atomic  ion-neutral  association.  This  reaction  involves  three  atoms, 
one  ionired  and  two  neutral,  as  shown  in  reaction  67. 
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X + X + X-X.  + x 


(67) 

Three  reactions  of  this  type  are  included  in  this  model  of  the  HgCl 
laser.  These  reactions  involve  the  formation  of  Ar+  and  Xe*  , 
and  are  steps  in  the  charge-transfer  process  that  creates  Hg+ 

The  first  theory  presented  is  Mahan's  theory  (Ref  26),  and  the  second 
theory  is  Smirnov's  theory  (Ref  27). 

Mahan's  theory.  Mahan's  theory  (Ref  26)  is  based  on  the 
assumption  that  the  ion  undergoes  a resonant  charge  transfer  with  a 
second  identical  neutral  particle  while  a third  arbitrary  particle  is 
at  some  distance  from  the  first  two  particles.  In  this  charge  transfer, 
a negligible  amount  of  kinetic  energy  is  lost.  The  second .particle 
that  is  now  charged  is  at  some  distance  from  a third  arbitrary 
uncharged  particle.  If  the  relative  kinetic  energy  of  the  third 
particle  and  the  newly  charged  second  particle  is  less  than  the  negative 
potential  energy  of  the  two  particles,  a bound  diatomic  molecule  is 
formed.  Since  no  kinetic  energy  is  transferred,  this  process  is 
called  a potential-energy  transfer  (Ref  26). 

Mahan's  next  assumption  is  that  the  potential  energy  between  the 
second  and  third  particles  is  the  ion-induced  dipole  type  (Ref  26). 
Hirschfelder  discusses  the  energy  of  ion-induced  dipole  type  reaction 
in  Reference  28.  Mahan  then  assumes  the  ion-induced  dipole  potential 
does  not  affect  the  second  and  third  particles  until  after  the 
resonant  charge  transfer  process  takes  place.  This  assumption  means 
the  charge  transfer  must  take  place  in  a short  time.  A fast  charge 
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transfer  occurs  if  the  first  and  second  particles  are  at  a large 
separation  distance  (Ref  26),  because  the  second  particle  spends  less 
time  in  the  sphere  of  interaction  of  particle  1 as  the  separation 
distance  of  the  particles  increases. 

These  assumptions  lead  to  a bound-molecular  ion  being  formed 
if  the  ion-induced  dipole  potential  energy  is  greater  than  the 
neutral  kinetic  energy  of  the  second  and  third  particles.  This 

•y 

• Ct  6^ 

condition  is  given  in  equation  68  where  (■^py)  is  the  ion-induced 

1 , 

dipole  potential  energy  and  (—  Mrv£)  is  the  relative  kinetic  energy 
of  the  second  and  third  particles. 


t M v* 
2 r r 


2r 


(68) 


In  equation  68,  M is  the  reduced  mass,  v the  relative  speed,  and 
r is  the  separation  distance  of  the  second  and  third  particles.  Also 
in  equation  68,  a is  the  polarizability  of  the  third  particle  and 
e is  the  electron  charge.  Solving  equation  68  for  r shows  that  a 
bound  molecule  is  not  formed  by  the  second  and  third  particle  unless 
equation  69  is  true. 

r = )/[*  (69) 

M v: 

r r 
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The  critical  radius , r%  , is  the  maximum  separation  distance  for  the 
second  and  third  particles  to  be  apart  and  still  form  a bound  molecule 
after  the  charge  transfer  takes  place. 

Mahan's  reaction-rate  constant  for  ion-neutral  association  is 
given  by  the  product  of  the  reaction-rate  constant  for  charge  transfer 
and  the  probability  that  the  third  particle  is  within  a sphere  of 
radius  r^  that  is  centered  at  the  second  particle  (Ref  26).  This 
condition  is  shown  in  equation  70  where  f(y  ) is  the  relative  speed 
distribution  and  is  the  reaction-rate  constant  for  charge  transfer 

of  the  first  and  second  particles. 


“l-N  • “CT  t L *«’<V  f(vr>  dV  (70) 
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Mahan  states  that  the  assumption  of  resonant  charge  transfer 
between  the  first  and  second  particle  is  not  strictly  true  if  the 
third  particle  is  close  enough  to  the  first  and  second  particles  to 
interfere  with  the  charge  transfer  reaction.  Mahan  also  states  that 
for  " typical  values  of  ionization  energy"  this  assumption  is  not 
incorrect  enough  to  affect  the  results  of  his  theory  (Ref  26) . 

If  the  speed  distribution  of  the  reacting  particle  is  assumed  to 
be  a Maxwell- Boltzmann  distribution  and  equation  69  is  substituted 
into  equation  “0,  the  reaction-rate  constant  for  ion-neutral  association 
is  given  by  equation  "1  and  71a  (Ref  26). 
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Mahan's  calculated  data  is  computed  by  multiplying  the  charge- 
transfer  cross-sections  of  Dalgarno  (Ref  29)  by  the  root-mean-square 
relative  speed.  Mahan  uses  the  polarizabilities  of  Pitzer  (Ref  30) 
for  the  calculated  values  (Ref  26).  Mahan's  method  of  calculation 
appears  to  be  reasonable,  because  his  calculated  values  are  all  within 
3 factor  of  three  and  most  are  within  a factor  of  two  of  the  experi- 
mental values  (Ref  26) . 

Using  Mahan's  method  of  calculation  and  the  data  in  references 
29  and  30,  a reaction-rate  constant  of  2.36  x 10’31  cm°/sec  is 
calculated  for  reaction  5 at  300°  K. 

Ar+  + 2Ar  -*•  Art  + Ar  (5) 

This  value  compares  necely  to  the  experimental  value  obtained  by 
Smith,  Dean,  and  Plumb  of  (2.30  : 0.10)  x 10" 31  cm5/sec  (Ref  31). 

The  value  that  is  calculated  for  reaction  5 at  533°  K is 

2.59  x 10  31  cmVsec  . The  value  that  is  calculated  for  reaction  7 

at  500°  K assuming  M is  Xe  is  5.3  x 10  - 1 cm°/sec 

Xe+  + Xe  + M - Xet  + M (M  = Xe,Ar)  (7) 

This  value  is  within  a factor  of  two  of  the  experimental  value  of 
(3.57  ± 0.17)  x 10  31  cm°/sec  that  is  reported  by  Smith,  Dean,  and 
Plumb  (Ref  31),  and  it  is  within  a factor  of  three  of  the  value 
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for  reaction  ' reported  by  Vitols  and 


(2.0  * 0.2)  \ lo"31  cm6 /sec 
Oskam  (Ref  32). 


The  reaction-rate  constant  for  reaction  ~ using  Ar  as  a third 
body  is  not  computable  using  Mahan's  theory  (Ref  2o)  and  is  an 
example  of  the  limited  applicability  of  Mahan's  theory.  This  reaction- 
rare  constant  cannot  be  calculated,  because  Mahan's  theory  assumes 
charge  transfer  between  the  first  and  second  particles  and  molecular 
formation  between  the  second  and  third  particle.  If  this  method 
of  formation  is  applied  to  reaction  ",  this  reaction  is  now  reaction  "2. 


Xe+  + Xe  + Ar  ■*  ArXe+  + Xe  (721 

The  products  in  reaction  72  can  be  formed  in  two  wavs,  first,  the 

+ 

charge  transfer  is  between  Xe  and  Xe  and  the  molecule  \rXe 

t 4 4* 

is  formed.  Second,  the  charge  transfer  is  between  Xe  and  Ar  and 
the  molecule  ArXe+  is  formed . Hie  second  method  is  not  included  in 
Mahan's  theory,  because  Mahan  requires  charge  transfer  between 
identical  particles. 

Smirnov 1 s theory . Smirnov's  theory  (Ref  2")  is  the  solution  to 
the  problems  presented  by  Malian's  theory.  Smirnov's  theory  presents 
methods  for  the  solution  of  reactions  73a,  b,  c,  and  d whore  X and 
V are  different  atoms. 
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(73a) 


X + X + V X j + ) 


l "3b) 
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X*  + X ♦ V ■+  XY+  +■  X 


("3c) 


X + X + X - X->  * X 


(73d) 


Smirnov  first  investigates  reaction  T'd. 

Smirnov's  theory  is  based  on  the  assumption  that  the  energy  of 
the  gas  particles,  which  is  proportional  to  the  temperature  of  the 
gas,  is  much  lower  than  the  dissociation  energy  of  a molecular  ion 
(Ref  27).  The  rate  of  change  of  atomic  ions  into  molecular  ions  is 
given  in  equation  "4  where  + is  the  molecular  ion  density, 

is  the  atomic  ion  density,  N'^  is  the  neutral  atom  density, 
and  Rj  ^ is  the  reaction-rate  constant  for  ion-neutral  association. 
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Smirnov  states  that  the  potential  energy  of  interaction  between  the 
atomic  ion  and  neutral  atom  is  determined  by  induced-dipole  polariration 
interaction  when  the  binding  energy  of  the  molecular  ion  is  on  the 
order  of  the  gas  temperature  (Ref  27). 

The  (Rj  ^ N^)  term  from  equation  74  is  examined.  This  term  is 
found  to  have  the  relationship  shown  in  equation  "3  (Ref  27). 


RI-N  NA  ’ NA<W>W 


The  (N ^ <cv>)  term  gives  the  frequency  that  a neutral  atom  creates 
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the  potential  to  be  bound  with  an  atomic  ion  to  form  a molecular  ion. 

The  W term  is  the  ratio  of  the  probability  that  the  molecular  ion 

gives  up  energy  in  a collision  with  a third  neutral  atom  and  creates 

a more  bound  molecular  ion  to  the  probability  that  the  molecular  ion 

gains  energy  on  the  order  of  the  temperature  in  a collision  with 

third  neutral  atom  and  dissociates.  As  a result  of  the  detailed 

balancing  principle  and  for  the  induced-dipole  polarization  inter- 
ae2  3/ 

action,  W * (-y)  4 where  a is  the  polarizability  of  the 

-,e2  1 A 

atom.  Smirnov  states  that  o is  proportional  to  (-y)  2 for 

the  polarization  interaction  cross  section  for  transition  of  the  ion 
into  weakly  excited  states  (Ref  27).  Equation  75  is  now  equation  76. 
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Na  <v> 


T 1 A 

Because  the  velocity  is  proportional  to  (— ) - for  a Maxwell-Boltzmann 

distribution,  equation  76  transforms  to  equation  77  where  M is  the 
mass  of  the  gas  atoms. 


,cte2  *A  Ae2'i3/4 


constant  (-rp)  " (-=-) 


Equation  77  is  true  for  reaction  7d. 

Smirnov  forces  equation  77  to  equate  to  1.1  x 10  31  cm°/sec  for 
reaction  73d  where  X is  He  and  the  temperature  is  300°  K.  This 
process  results  in  a value  for  the  constant  of  approximately  42  in 
equation  77  (Ref  27) . 
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The  results  using  Smirnov's  theory  (Ref  27)  are  not  as  good  as 
those  using  Mahan's  theory  (Ref  26).  Using  the  data  from  Titter 
(Ref  50),  a reaction-rate  constant  of  5.24  x 10_il  cm°/sec  is 
calculated  for  Ar  ion-neutral  association  in  reaction  5 at  a 
temperature  of  500°  K.  This  value  is  more  than  a factor  of  two  greater 
than  both  the  value  of  2.56  x 10  31  cm°/sec  from  Mahan's  theory  and 
the  experimental  value  of  (2.5  ± 0.10)  x 10  31  cm6/sec  (Ref  51)  at 
the  same  500°  K temperature.  The  value  that  is  calculated  for  reaction  7 
at  500°  K with  Xe  as  the  third  body  using  Smirnov's  theory  is 
3.38  x 10  cm" /sec  . This  value  is  more  than  a factor  of  two  greater 
than  the  experimental  value  of  (5.57  ± 0.17)  x 10  31  cm6/sec  (Ref  51). 
Mahan's  theory  gives  a value  within  a factor  of  two  of  the  experimental 
value.  Smirnov's  theory  does  give  results  within  a factor  of  three 
of  the  experimental  results  for  reactions  5 and  7. 

Smirnov  (Ref  27)  modifies  equation  ”7  using  Thomson ' s theory 
(Ref  55)  to  give  reaction-rate  constants  valid  for  reaction  "5a,  b, 
and  c.  The  assumption  necessary  for  this  modification  is  that  the  mass, 
m , of  particle  Y is  less  than  the  mass,  'I  , of  particle  X 
The  reaction-rate  constants  for  reactions  ~3a,  b,  and  c are  given  in 
equations  77a,  b,  and  c,  respectively  (Ref  27). 
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Smirnov  (Ref  27)  notes  the  value,  1.5  x 10~31  cm°/sec,  that  is 
calculated  for  reaction  78  at  approximately  400°  C is  close  to  Biondi’s 
experimental  value  (Ref  54)  of  1.7  x 10  31  cm6/sec. 


Hg+  + Hg  + He  Hg2  + He 


(78) 


This  close  agreement  for  one  reaction  does  not  imply  that  a greater 
accuracy  than  the  factor  of  three  that  is  applied  to  equation  77  is 
applied  to  equations  77a,  b,  and  c. 

The  value  that  is  calculated  for  reaction  7 using  Ar  as  a third 
body  and  Smirnov's  equation  77b  is  4.96  x 10  31  cm°/sec  3t  535°  K. 

The  value  that  is  calculated  for  reaction  7 using  Xe  as  a third 
body  and  Smirnov's  equation  77  is  5.77  x 10  31  cm5/sec  at  533°  K. 
Reaction  7a  is  reaction  7 with  Ar  as  the  third  body,  and  reaction 
7b  is  reaction  7 with  Xe  as  the  third  body. 

Xe+  + Xe  + Ar  - Xe?  + Ar  (7a) 

Xe+  + Xe  + Xe  -*■  Xe?  + Xe  (7b) 


Conclusion . Mahan's  theory  (Ref  26)  appears  to  predict  reaction- 
rate  constants  to  within  a factor  of  two  of  known  experimental  results. 
This  theory  is  limited  in  applicability,  because  Mahan  requires  a 
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specific  order  of  charge  transfer  and  interaction.  Smirnov's  theory 
(Ref  27)  applies  to  more  types  of  reactions  than  Mahan's  theory.  The 
values  predicted  by  Smirnov's  theory  are  probably  within  a factor  of 
three  of  known  experimental  values.  Therefore,  Mahan's  theory  gives 
more  accurate  results  for  a limited  number  of  reactions,  and  Smirnov's 
theory  gives  results  for  some  types  of  ion-neutral  association  reactions 
not  computable  from  Mahan's  theory. 

The  following  values  are  used  in  the  model  for  ion-neutral 
reaction-rate  constants.  The  experimental  value  (Ref  31)  of 
2.5  x 10~31  cm5/sec  is  used  in  the  model  for  reaction  5.  The  experi- 
mental value  (Ref  51)  of  3.6  x 10  31  cm5/sec  is  used  for  reaction  7b. 

Reaction  7b  is  reaction  7 with  Xe  as  the  third  bodv. 

* 

I 

Xe+  + Xe  + Xe  -*■  Xe*  + Xe  (7b) 

The  value  that  is  calculated  using  Smirnov's  equation  77a, 

4.96  x 10~31  cm3/sec  , is  used  for  reaction  7a.  Reaction  7a  is 
reaction  7 with  Ar  as  the  third  body. 

Xe*  + Xe  + Ar  ■+  Xe?  + Ar  (7a) 

The  next  section  is  a summary  of  the  charge  transfer  reactions 
in  this  model. 


Charge  Transfer  Reactions 

Introduction.  The  charge  transfer  processes  that  are  included 
in  this  section  are  the  near-symmetric  charge  transfers  in  reactions 
6 and  8 of  this  model.  The  energy  of  the  particles  is  approximately 
0.07  eV.  The  charge  transfer  between  Xe?  and  Hg  in  reaction  8 
is  a near-symmetric  charge  transfer  (Ref  35) . The  general  form  of 
reactions  to  be  considered  in  this  section  is  shown  in  equation  79 
where  X is  an  atom  and  Y is  a different  atom. 

*2  + Y - Y+  + 2X  (79) 

Two  theories  are  presented.  The  first  theory  is  that  of 
Gioumousis  and  Stevenson  (Ref  36).  The  values  generated  by  using 
Gioumousis  and  Stevenson's  theory  are  in  excellent  agreement  with 
the  experimental  values  reported  by  Bohme  et  al.  (Ref  57).  The 
second  theory  is  a brief  presentation  of  Rapp  and  Francis'  theory 
for  the  low-velocity  region  (Ref  38) . The  particles  in  reactions 
6 and  8 are  in  the  Rapp  and  Francis  low-velocity  region  that  is 
determined  by  the  approximate  upper  limit  v : (105/Mr  ' 2)  cm/sec 
where  M is  the  reduced  mass  of  the  reacting  particles  in  amu.  Since 
the  velocity  of  the  gas  particles  is  proportional  to  the  square  roots 
of  the  temperature  of  the  gas,  the  Rapp-Francis  low  velocity  region 
is  defined  by  the  gas  temperature  and  the  reduced  mass  of  the  reacting 
particles. 

Gioumousis  and  Stevenson's  theory.  Gioumousis  and  Stevenson's 
theory  is  based  on  the  fact  that  the  long-range  potential  between  an 
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ion  and  a neutral  particle  is  given  by  equation  SO  which  is  the  ion- 
mduced  dipole  polarization  potential  (Ref  56). 


U » 


(SO) 


In  equation  SO,  U is  the  potential  energy,  e is  the  electronic 
charge,  is  the  polarizability  of  the  neutral  particle,  and  r is 

the  separation  distance  of  the  ion  and  the  neutral  particle.  Gioumousis 
and  Stevenson  use  Langevin's  (Ref  59)  equation  below  describing  the 
orbits  for  the  ion-induced  dipole  potential  (Ref  56).  Equation  SI 
describes  these  orbits  where  b is  the  impact  parameter. 
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In  equation  SI,  M is  the  reduced  mass  of  the  reacting  particles  and 

v is  their  relative  velocity.  A typical  set  of  trajectories  is 
r 

included  in  Figure  S.  The  important  facts  in  Figure  S are  that  for 

b < b^  the  orbits  pass  through  the  center  point  and  for  b t b^  the 

orbiting  particles  are  never  closer  than  (b^/ZT) 

The  next  step  in  the  development  of  this  theory  is  the  supposition 

that  a critical  radius.  r„  , for  a given  reaction  to  take  place 

exists  (Ref  56).  '11115  critical  radius  is  such  that  the  probability 

of  the  reaction  taking  place  is  zero  if  the  particle  passes  outside 

the  critical  radius  and  one  if  the  particle  passes  inside  the  critical 

radius.  Therefore,  if  r d b AT  the  cross  section  of  interaction 

c o 

is  given  in  equation  S2. 


Fig.  S.  A typical  family  of  trajectories  as  a function  of  the  impact 
parameter  b.  The  dotted  trajectory  is  the  critical  one  for  b = b 
and  approaches  the  circle  r = b^/ZlT  . Only  trajectories  with 

b < b will  enter  the  reaction  sphere  if  it  has  radius  r less  than 
0 c 

this,  i.e.,  r < b /»'?.  Thus,  the  cross  section  for  entrance  into 
c o 

the  reaction  sphere  depends  on  b but  not,  within  limits,  r . On 

0 c 

the  assumption  that  all  molecules  which  enter  the  reaction  sphere 
do  react,  the  same  may  be  said  of  the  interaction  cross  section  (.Ref  ?6j 


The  larger  cross  sections  are  given  for  interactions  between 
particles  that  maximice  b^  . This  maximization  is  done  by  particles 
with  3 large  polarizability,  small  reduced  mass,  and  small  velocities. 
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Since  velocities  are  often  described  in  terras  of  a velocity  distribution, 
some  higher  velocities  are  certain  to  exist  so  that  the  condition 
r < b^/*?  is  not  true.  If  the  number  of  particles  exceeding  this 
criteria  is  kept  to  3 minimum,  this  theory  is  still  relatively  valid 
(Ref  56). 

The  cross  section  in  equation  82  multiplied  by  the  relative 
velocity  of  the  two  particles  is  a constant  shown  in  equation  83. 
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Therefore,  the  expectation  value  of  <v  a(v  )>  is  the  constant  given 
in  equation  85  for  ions  and  neutral  particles  in  a Maxwel 1-Boltrmann 
distribution.  This  constant  is  the  reaction-rate  constant  for 
symmetric  charge  transfer. 

Rapp  and  Francis'  theory.  The  Rapp  and  Francis  theory  (Ref  3S) 

1 0 5 

for  the  low-velocity  region,  v f ( — r.)  , is  a modification  of  the 
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Gioumousis  and  Stevenson  theory  (Ref  36).  Rapp  and  Francis  say  that 
the  charge  exchange  probability  is  0.5  for  symmetric-charge-transfer 
reaction  resulting  from  impact  collisions  (Ref  58).  By  multiplying 
this  factor  of  0.5  times  the  Gioumousis  and  Stevenson  cross  section 
of  reaction  given  in  equation  82,  Rapp  and  Francis  give  a lower  bound 
cross  section  for  symmetric  charge  transfer  shown  in  equation  84. 
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Using  the  approximate  value  of  10  -1*  cm3  for  the  polarizability  of 

particle  B allows  Rapp  and  Francis  to  approximate  the  cross  section. 

This  approximate  cross  section  is  given  in  equation  85  where  M is 

r 

in  amu  and  v is  in  cm/sec. 
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Experimental  data.  The  charge-transfer  experiments  of  Bohme 
et  al.  (Ref  37)  provide  excellent  results  for  comparison  to  the 
theoretical  reaction-rate  constants  that  are  calculated  using 
Gioumousis  and  Stevenson's  theory  (Ref  36).  For  example,  the 
theoretical  value  for  reaction  86  is  5.73  x 10  cmJ/sec  at 

200°  K,  and  the  experimental  value  is  7.5  x 10~13  cm3/sec  (Ref  57). 
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Kr+  + 
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(86) 


Most  of  Bohme 's  data  (Ref  37)  is  within  a factor  of  two  of  the 
theoretical  values.  The  values  not  within  a factor  of  two  of  the 
theoretical  values  are  less  than  the  theoretical  values  and  are  best 
explained  as  reactions  that  are  dominated  by  asymmetric  charge  transfers 
(Ref  37) . The  experimental  data  of  Maier  for  particles  with  greater 
than  thermal  energies  does  not  agree  with  Bohme's  data  (Ref  401. 
Therefore,  Gioumousis  3nd  Stevenson's  theory  can  be  used  for  calculating 
theoretical  values  of  charge-transfer  reaction-rate  constants  for 
particles  with  thermal  energies.  These  theoretical  values  can  be 
considered  to  be  within  a factor  of  two  of  the  upperlimit  reaction- 
rate  constants.  If  the  charge  transfer  is  asymmetrical,  the 
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theoretical  value  places  an  upper  limit  on  the  reaction-rate 
constant . 

The  theoretical  value  for  reaction  6 in  this  model,  Ari 

transferring  charge  to  Xe  , is  6.65  x 10”  cm3/sec  . The 

theoretical  value  for  reaction  S in  this  model,  Xe;  transferring 

charge  to  Hg  , is  3.11  x 10  cnr/sec  • The  polar inability 

of  Xe  from  Fitter  (Ref  30)  is  used  in  the  calculation  for  reaction  6. 

-24  7 

. An  estimated  value  of  2 x 10  cnr  is  used  for  the  polarizability 

of  Hg  in  reaction  3. 

Bohme  et  al.  (Ref  57)  also  use  Gioumousis  and  Stevenson's 
theory  (Ref  36)  to  calculate  theoretical  values  for  reactions 
similar  to  reaction  37. 

ArKr+  + Kr  - Kr£  + Ar  (87) 

The  general  form  of  reaction  S7  is  shown  in  reaction  88  where  X and 
Y are  different  atoms. 

XY+  ♦ Y -v  Yj  + X (SS) 


This  theoretical  value  for  reaction  S7  is  5.21  x 10  10  cmVsec  » 
and  the  experimental  reaction-rate  constant  is  3.2  x 10  ‘ cnr/sec 
(Ref  37).  The  theoretical  values  that  are  calculated  tor  this  type 
reaction  are  generally  high  by  a factor  of  three  or  less  (Ref  57). 
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Conclusion.  The  Gioumousis  and  Stevenson  theory  (Ref  36) 
results  in  reaction-rate  constants  within  a factor  of  two  for  reactions 
similar  to  reaction  "9  for  particles  with  thermal  energies.  If  the 
reaction  is  an  asymmetric  charge  transfer,  the  theoretical  value  is 
an  upper- limit  reaction-rate  constant.  The  Rapp-Francis  theory 
(Ref  38")  at  low-velocities  is  not  as  accurate  as  the  Gioumousis  and 
Stevenson  theory.  The  Rapp-Francis  criteria  for  determining  low- 
velocity  reactions  is  probably  a reasonable  criteria  for  determining 
the  applicability  of  the  Gioumousis  and  Stevenson  theory.  The 
fact  that  the  Gioumousis  and  Stevenson  theory  Joes  not  predict  accurate 
reaction-rate  constants  for  energies  greater  than  thermal  energies 
is  pointed  out  in  a comparison  of  Bohme's  thermal  energy  data  (Ref  37) 
to  Maier's  data  (Ref  40)  for  energies  greater  than  thermal.  Therefore, 
the  Gioumousis  and  Stevenson  theory  gives  excellent  reaction-rate 
constants  for  the  charge-transfer  reactions  that  are  similar  to 
reactions  6 and  3 with  thermal  energy  reactants. 

The  next  section  includes  discussion  about  the  three  ionic 
recombination  processes  in  this  model. 

Ion- Ion  Recombination 

Introduction.  Ion-ion  recombination  reaction  is  the  process  by 
which  a singly  ioniced  ion  combines  with  a negative  ion  that  has 
one  extra  electron.  The  four  theories  investigated  in  this  section 
are  J.  J.  Thomson's  theory  (Ref  53),  Langevin's  theory  (Ref  41), 
Natanson's  theory  (Ref  41),  and  Wadehra  and  Rardsley's  theory  (Ref  42). 
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Thomson's  theory.  Thomson's  theory  of  ionic  recombination 
(Ref  33)  is  based  on  two  assumptions.  First,  all  gases  are  in  thermal 
equilibrium.  Second,  if  two  oppositely  charged  ions  pass  within  a 
certain  distance  of  each  other,  called  the  trapping  radius,  these  ions 
recombine  (Ref  45).  This  development  is  similar  to  that  in  Reference 
43. 

The  assumption  that  all  ions  are  in  thermal  equilibrium  with  the 
gas  means  the  average  kinetic  energy  of  the  ions  is  equal  to  the 
average  energy  of  any  particle  in  the  gas.  This  relationship  is 
shown  in  equation  89  where  mr  are  the  ionic  masses,  v~  are  the 
root  mean  square  velocities  of  the  ions,  T is  the  temperature  of 
the  gas  in  °k  , and  k is  the  Boltzmann  constant. 


1 + +2 
— m v 
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If  the  distance  between  the  oppositely  charged  ions  is  large, 
then  equation  90  is  true  where  M is  the  reduced  mass  of  m and 
m+  , and  vnD  is  the  root-mean-square  velocity  of  the  relative 
velocities  of  the  ions. 
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The  next  assumption  is  that  the  potential  energy  between  the  two 
ions  is  only  Coulombic.  This  assumption  means  that  as  the  ions  are 
separated  by  greater  and  greater  distances,  the  Coulomb  potential 
energy,  -e2/r  , approaches  zero. 
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rhomson's  assumption  that  loads  to  the  idea  of  a trapping 
radius  is  that  two  ions  recomhine  if  their  potential  energy  and 
kinetic  energy  sum  to  a negative  number.  Equation  $0  shows  the 
average  kinetic  energy  of  the  two-ion  system  is  ’ kT  , which  is 

positive.  Since  the  two  ions  gain  kinetic  energy  equal  to  their  loss 
of  potential  energy  as  the  two  ions  approach  each  other,  the  only  way 
the  kinetic  energy  is  ever  less  than  the  magnitude  of  the  potential 
energy  is  for  kinetic  energy  to  he  given  up  in  a collision  with  a 
third  body.  The  kinetic  energy  of  the  ion  starts  at  ’ kT  and 
increases  as  it  approaches  the  other  ion.  If  an  ion  is  returned  to 
an  equilibrium  temperature  by  a collision,  then  that  ion  is  trapped 
when  c'/r  1 2 ^T  • Tit®  distance,  r , at  which  this  occurs  is 

ca 1 1 ed  t he  t rapp i ng  rad ius . 


Thomson's  next  assumption  is  that  each  ion  is  surrounded  by  a 
sphere  with  a radius  equal  to  the  trapping  radius.  The  number  of 
positive  ions  entering  a sphere  of  interaction  around  a negative  ion 

, +-1  -;!/■>  5 

with  radius  r_  is  t r;,.n  (v  ~ + v ) - , where  nr*  is  the  circle 

T F + m m 1 

of  interaction,  n is  the  number  density  of  positive  ions,  and 

+:  .2  i/,  „ 

(v  * v ) ' - is  the  average  relative  velocity  ot  the  ions.  The 
mm 

number  of  negative  ions  entering  a positive  sphere  is 
, *2  -21/ 

ir_;n  iv  + v 1 '2  where  n Is  the  number  densitv  of  negative 
r - m m 

ions,  hot  w*  be  the  probabilities  that  a positive  or  negative  ion 
entering  an  appropriate  sphere  of  interaction  collides  with  another 
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particle  and  is  thermaliced  in  that  sphere.  The  number  of  recombina- 
tions per  cubic  centimeter  per  second  is  then  given  in  equation  92. 
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An  approximate  value  for  w“  can  be  calculated  using  classical 

mechanics  and  ignoring  the  charges  of  the  two  ions.  First,  the 

probability  that  a particle  travels  a distance  x in  a gas  and  does 

xA 

not  collide  with  another  particle  is  e”  , where  \ is  the  mean 
free  path  of  the  moving  particle.  The  mean  free  path  is  a function 
of  the  total  gas  density  and  is  not  a function  of  only  n+  and  n_  . 
Second,  the  distance  x traveled  by  a particle  in  the  sphere  is  given 
by  equation  93. 


x = 2rT  cos  9 


(931 


The  angle  9 is  formed  by  the  intersection  of  the  trajectory  line 
and  the  line  formed  by  the  point  of  entry  to  the  circle  and  the  center 
of  the  circle.  Third,  the  probability  that  a particle  enters  the  circle 
of  interaction  between  9 and  3+d9  is  sin  9 (Ref  44).  Therefore, 
the  probability  that  a particle  enters  a circle  with  radius  r^  and 
collides  with  another  particle  is  given  in  equation  94. 


Probability  of  ■> 
Collision 


1 - J e“rTCOse/xsin  0 d 9 
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The  actual  probabilities  of  collision  are  greater  than  those 
calculated  in  equation  94,  because  the  effect  of  the  difference  in 
charge  is  ignored. 

+ 

Substituting  w"  in  equation  92  and  dividing  by  (n+n_) 
results  in  the  reaction-rate  constant  for  ion-ion  recombination  given 
in  equation  95  where  \+  are  the  mean  free  paths  of  the  positive  and 
negative  ions. 
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If  the  ions  and  the  gas  are  assumed  to  be  in  a Maxwell-Boltomann 

+2  .2  1 / qi.t  1 / - 

distribution,  (v  + v 1 1 equals  (-£— ) - 

m m 1 

r 

Two  cases  can  be  derived  for  equation  95.  The  first  case  is  for 

_ o r / \ 

low  pressure  gases.  At  low  pressures,  \ are  large  and  e “ T ± 

can  be  estimated  by  the  first  three  terms  of  the  series  representation 

— y + 

for  e 7 where  y is  small.  The  term  w is  approximately  ry/^+ 

and  w is  approximately  r^,/\  . The  reaction-rate  constant  for  low 

pressure  ion- ion  recombination  is  given  in  equation  96. 
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The  second  case  is  for  high  pressure  gases.  At  high  pressures, 

_2r  /X 

are  small  and  e " T ± approaches  one.  The  reaction-rate  constant 
for  high  pressure  is  given  in  equation  97. 
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Substituting  equation  91  into  equations  96  and  97  results  in 
equations  96a  and  97a. 
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Equations  96a  and  97a  contain  some  interesting  information. 

The  low  pressure  equation  96a  indicates  that  at  a constant  temperature 
the  reaction-rate  constant  increases  with  pressure.  The  high  pressure 
equation  97a  shows  that  the  ion-ion  recombination  reaction-rate 
constant  is  independent  of  pressure  at  constant  temperature. 

Several  other  comments  are  made  in  Reference  45  concerning  the  low 
pressure  equation  97’a,  the  main  point  being  th3t  the  Thomson  low 
pressure  reaction-rate  constant  for  ion-ion  recombination  is  a good 
approximation  for  pressures  less  than  1 atmosphere. 

The  fact  that  the  Thomson  high  pressure  reaction-rate  constant 
is  independent  of  pressure  is  in  disagreement  with  experimental  data 
(Ref  45).  The  reason  for  Thomson's  generally  high  prediction  of  the 
recombination  constant  is  that  no  consideration  is  given  to  the 
possibility  of  multiple  collisions  inside  the  sphere  of  interaction 
knocking  a possible  reacting  particle  out  of  the  sphere  of  interaction. 


' 
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Also,  no  consideration  is  given  to  dissociative  collisions  between 
molecules  (Ref  46).  These  facts  cause  the  Thomson  theory  generally 
to  overestimate  the  high  pressure  reaction-rate  constant  for  ion- ion 
recombination. 


Langevin's  theory.  Langevin's  theory  is  valid  only  at  very  high 
pressures  and  is  based  on  the  assumption  that  the  attractive  forces 
between  two  neighboring  particles  is  Coulombic.  The  positive  and 
negative  ions  drift  toward  each  other  at  a rate  determined  by  their 
mobilities  K~  . The  relative  drift  velocity,  v^  , at  a separation 
distance  r of  the  two  ions  is  ( K+  + K~)  e2/r2  , where  the  effective 

electric  field  on  the  two  ions  is  e2/r2  . If  each  negative  particle 

is  enclosed  in  a sphere  of  radius  r , the  number  of  positive  particles 
entering  the  sphere  is  given  by  the  product  of  the  drift  velocity, 
the  surface  area  of  the  sphere,  and  the  number  density  of  negative 
ions.  The  number  of  recombinations  is  shown  in  equation  9S. 


• — + - 

n , . = 4re-  (K  + < ) n n 

recombinations  v 


gevin's  reaction- rate  constant  for  ion- ion  recombination  is  given 
in  'equation  99. 


(98) 


Rj.il  = (K+  + <’) 


(99) 


Langevin's  theory  does  agree  with  Machler's  data  at  5 to  12 
atmospheres  of  pressure  (Ref  41).  Langevin's  theory  also  predicts 
a decrease  in  the  reaction-rate  constant  as  pressure  increases. 
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The  major  problem  that  is  found  for  using  Langevin's  theory  for 
quick  calculations  for  ion-ion  recombination  is  the  lack  of  data 
for  r . 

Natanson's  theory.  Natanson's  theory  (Ref  471  fills  the  gap 
between  high  and  low  pressures.  Unfortunately,  Natanson's  theory 
requires  more  data  than  either  Thomson  or  Langevin  to  calculate 
the  reaction-rate  constant  for  ion-ion  recombination.  The  require- 
ment by  both  Langevin  and  Natanson  for  data  that  is  difficult  to 
obtain  in  calculating  the  reaction-rate  constants,  forces  the  use 
of  Thomson's  theory  for  approximate  calculations  of  ion-ion 
recombination  reaction-rate  constants. 

Wadehra  and  Bards  ley's  theory.  IVadehra  and  Bards  ley's  theory 
(Ref  42)  is  an  extension  of  Natanson's  theory  to  allow  the  masses 
and  collision  properties  of  the  positive  and  negative  ions  to  be 
different.  By  allowing  for  different  collision  properties  for  the 
positive  and  negative  ions,  a trapping  radius  is  developed  for 
the  positive  ion  that  is  different  from  the  trapping  radius  of  the 
negative  ion.  Wadehra  and  Bards  ley  cite  several  sources  of  data 
for  their  calculations. 

Conclusion.  The  use  of  Thomson's  theory  may  not  prove  to  be  as 
inaccurate  as  expected  in  the  calculations  for  ion-ion  recombination 
in  the  HgCl  laser.  First,  the  ion-ion  recombinations  in  the  HgCl 
laser  are  atom-atom  recombinations  and  not  molecule-molecule 
recombinations.  This  fact  eliminates  the  problem  of  dissociative 
collisions  considered  by  Brueckner  (Ref  46).  Second,  the  operating 
pressure  forth?  HgCl  laser  being  investigated  in  this  model  is 
approximately  3.17  atmospheres.  This  pressure  is  not  the  extremely 


high  pressure  at  which  Langevin's  theory  is  accurate.  This  pressure 
is  probably  such  that  the  reaction-rate  constant  is  closer  to  its 
maximum  value.  Evidence  of  this  fact  is  given  in  a theoretical 
calculation  done  by  M.  R.  Flannery.  Flannery's  calculations  for 
Hg+  and  Cl  ionic  recombination  result  in  a maximum  recombination 
constant  of  2.4  x 10  6 cm3/sec  at  N/N^  of  1.6  . This  reaction-rate 

constant  drops  rapidly  to  1 x 10  ° cm3/sec  at  N/NL  = 9.0  (Refs  48,  49). 
In  the  term  Nt/N^  » N is  the  density  of  the  third  body  Ar  , 
and  >y  is  Loschmidt's  number  that  is  adjusted  for  temperature.  If 
the  assumption  is  made  that  the  sum  of  the  number  densities  of  Xe 
and  Ar  in  this  model  is  a good  approximation  for  N , an  ionic 

-6  5 

recombination  reaction-rate  constant  of  2.15  x 10  cnr/sec  is 
found  from  Flannery's  data  to  apply  to  this  model.  Flannery's  data 
is  included  in  Appendix  C. 

Wadehra  and  Bards  ley  (Ref  42)  predict  an  ionic  recombination 
rate  of  2.3  x 10'6  cm3/sec  for  Ar+  and  Cl”  . This  rate  is 
computed  for  a gas  temperature  of  300°  K and  a pressure  of  approximately 
3.1  atmospheres.  The  predicted  ionic  recombination  rate  for  Hg+ 
and  Cl”  at  the  same  temperature  and  pressure  is  1.8  x 10”°  cm3/sec 

Three  ion-ion  recombination  reaction-rate  constants  based  on 
Thomson's  theory  are  used  in  this  model  of  the  HgCl  laser.  These 
values  result  from  calculations  using  equation  97a  and  are  presented 
in  Table  1. 


_ _ 
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Table  1 

Ion-Ion  Recombination  Constants  for  the  HgCl  Laser 


Reaction-Rate  Constants 
Reaction (cm3/sec) 


Ar+ 

+ 

+ Cl* 

- ArCl* 

2.13  x 10"6 

-6 

« 

Xe 

+ 

+ Cl" 

-*•  XeCl* 

1.75  x 10 

Hg 

+ Cl* 

-*  HgCl* 

1.6S5  x 10* 

These  approximate  values  calculated  for  the  HgCl  laser  are 
probably  within  a factor  of  three  of  the  actual  values.  This 
hypothesis  of  close  agreement  is  based  on  the  excellent  agreement 
of  the  value  calculated  using  Thomson's  model  with  Flannery's  Hg+ 
and  Cl  ionic  recombination  (Ref  48)  at  similar  operating  pressures. 
IVadehra  and  Bardsley's  calculations  (Ref  42)  also  provide  excellent 
data  3nd  agreement  with  values  calculated  using  Thomson's  theory. 

The  next  section  includes  a discussion  of  the  harpooning  theory 
as  it  is  applied  to  this  model. 

Harpooning  and  Similar  Reactions 

Introduction.  The  harpooning  reaction  shown  in  reaction  100 
usually  involves  an  alkali  atom,  M , harpooning  a diatomic  halogen 
molecule,  X., 


M + X,  MX  + X 

The  theory  that  is  presented  in  this  section  is  based  on  Herschbach' 
(Ref  50)  review  of  the  harpooning  reaction.  The  experimental  data 


(100) 
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(Ref  6)  that  is  included  in  this  section  shows  that  excited  states  of 
some  atoms  are  similar  to  alkali  atoms  and  that  these  excited  states 
are  involved  in  reactions  similar  to  harpooning  reactions.  Reactions 
such  as  Ar*  and  Xe*  reacting  with  Cl2  to  form  ArCl*  and  XeCl* 
are  also  included  in  this  section. 

Theory.  The  basic  assumption  that  is  made  by  Magee  (Ref  51) 

is  that  the  potential  energy  curves  of  MX  and  (M  + X)  cross  at 

some  point  r^  that  is  called  the  capture  point.  This  crossing  is 

shown  in  Figure  9 for  KBr  (Ref  50).  At  the  internuclear  separation 

distance  r , the  Coulombic  attraction  forces  are  dominant.  The 
c 

capture  point,  r,  , can  be  determined  by  finding  the  energy  that 
is  required  to  form  the  ion  pair.  The  energy  that  is  required  for 
this  reaction  is  essentially  determined  by  the  amount  of  energy  that 
is  required  to  strip  an  electron  from  the  alkali  atom  minus  the  energy 
that  is  given  up  by  the  halogen  atom  after  the  halogen  atom  captures 
the  electron  of  the  alkali  atom.  Eequation  101  (Ref  50)  shows  the 
potential  energy  equal  to  the  energy  that  is  required  for  the  harpoon 
reaction. 


e2/rc  = I (M)  - E(X)  (101) 

In  equation  101,  e is  the  electronic  charge,  I (M)  is  the  ionitation 
potential  of  the  alkali  atom,  and  E(.X)  is  the  electron  affinity  of 
the  halogen  atom. 
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Fig.  9.  Potential  energy  curves  for  an  alkali  halide  molecule  (drawn 
for  KBr)  showing  the  "zeroth-order  crossing"  of  the  ionic  and  covalent 
states  (Ref  50) . 


Krause  et  al.  (Ref  6)  hypothesize  that  the  Hg*  (3P2)  reaction 
with  CI2  to  form  HgCl*  is  similar  to  the  harpoon  reaction.  A 
calculation  using  an  ionization  potential  of  4.974  eV  (Ref  6)  for 
Hg*  (3P2)  and  an  electron  affinity  of  2.19  eV  for  CI2  (Ref  21) 

O 

results  in  a reaction  cross  section  of  S4.5  A2  . This  cross  section 


is  calculated  by  using  irr2  as  the  reaction  cross  section.  This 


calculated  cross  section  is  in  good  agreement  with  the  experimental 
cross  section  for  this  reaction  that  is  determined  by  Kraus  et  al. 
of  (90  ± 25)  A2  (Ref  6) . 

Experimental  work.  By  multiplying  the  experimental  cross  section 


of  Kraus  et  al.  (Ref  6)  by  the  average  relative  velocity  of  the  Hg1 


10 


and  Cl2  in  this  model,  a reaction-rate  constant  of  4. IS  x 10  J cm3/sec 


I 


I 
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is  calculated  for  IlgCl*  formation  by  Hg*  harpooning  of  CI2  . These 
gas  particles  in  the  model  are  assumed  to  be  in  a Maxwell -Boltzmann 
distribution  at  a temperature  of  533°  K. 

Reactions  that  are  included  in  this  model  that  are  similar  to  the 
HgCl*  formation  by  Hg*  and  CI2  are  the  formation  of  ArCl*  by  Ar* 
and  Cl2  and  the  formation  of  XeCl*  by  Xe*  and  CI2  . The  reaction- 
rate  constants  for  these  reactions  are  calculated  assuming  these  gas 
particles  are  in  a Maxwell -Boltzmann  distribution  at  533°  K.  The  XeCl* 
formation  reaction-rate  constant  is  calculated  using  the  quenching  cross 
section  for  Xe*  by  Cl2  that  is  reported  by  Velazco  and  Setser  to  be 

O 

210  ,\~  (Ref  52) , Velazco  and  Setser  (Refs  52,  53)  report  that  almost 
100%  of  the  quenching  of  Xe*  by  Cl2  includes  the  formation  of  XeCl* 

Therefore,  the  total  quenching  cross  section  is  multiplied  by  the  average 
relative  velocity  to  get  a reaction-rate  constant  of  1.04  x 10’9  cm3/sec 
for  formation  of  XeCl*  by  Xe*  and  Cl2  reacting.  The  ArCl*  formation 

O 

reaction-rate  constant  is  calculated  by  using  the  95  A cross  section  that 
is  reported  by  Piper  et  al.  (Refs  54,  55)  to  be  the  Ar*  quenching 
cross  section  by  Cl2  . Golde  and  Thrush  (Ref  7)  state  that  more  than 
10%  of  Ar*  quenching  by  Cl2  results  in  ArCl*  formation.  Gundel 
et  al.  (Ref  56)  state  that  approximately  50%  of  Ar*  quenching  by  Cl2 
results  in  ArCl*  formation.  Therefore,  taking  50%  of  the  total  quenching 
cross  section  and  multiplying  that  number  by  the  average  relative  velocity 
of  Ar*  and  Cl2  results  in  a reaction-rate  constant  of  3.16  x 10" 10 
cm3/sec  for  formation  of  ArCl*  by  Ar*  reacting  with  Cl2 
Conclusion . The  harpoon  reaction  does  seem  to  explain  the 
formation  of  HgCl*  by  Hg*  reacting  with  Cl2  . Unfortunately, 
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this  determination  can  be  made  only  after  the  experimental  data 
is  gathered  for  the  HgCl*  function  cross  section. 

The  reamining  sections  are  all  brief  discussions  about  the 
development  of  several  reaction-rate  constants  that  are  based  on 
experimental  data. 


Secondarv-Electron  Excitation  of  Ar,  Xe,  and  Hs 


Introduction.  Secondarv-electron  excitation  of  an  atom  is  a 


reaction  of  the  form  in  reaction  102. 


e + Y -*■  Y * ♦ e 


(102) 


The  Y*  is  an  excited  state  which  is  long-lived  and  is  called  a 
metastable  state. 

The  reaction-rate  constant  for  excitation  of  Ar  , Xe  , and 
Hg  by  secondary  electron  impact  are  based  on  theoretical  and  experi- 
mental cross  sections.  The  cross  sections  are  approximated  by 
straight  lines  of  the  form  in  equation  103,  where  c is  electron  energy 
in  eV  and  a is  the  resultant  cross  section.  The  units  of  a and 
b are  cnr/eV  and  cm-  , respectively. 

o = ae  + b (103) 


The  expectation  value  of  <ov>  is  calculated  for  a Maxwell-Boltzmann 

energy  distribution  of  electrons  with  4 eV  average  energy.  Using 
“V  1 / 

v = (— ) 2 allows  the  integration  limits  to  be  zero  energy  to 

m 

infinite  energy  as  shown  in  equation  104. 


’1 


excitation 


The  integration  yields  equation  10S  which  must  he  evaluated  for  the 


ero  to  infinite  energy  limits 


Ar  and  Xe.  The  Ar  and  Xe  cross  sections  for  electron- impact 


excitation  to  metastable  states  are  the  experimental  values  obtained 


by  Schaper  and  Schiebner  (Ref  57) . The  Ar  cross  section  is 


approximated  by  equation  lOo 


6.6666  x 10 


The  Xe  cross  section  is  approximated  by  equation  107 


Substitution  of  equations  106  and  107  into  equation  105  and  evaluation 


over  the  proper  limits  yields  electron-impact  excitation  reaction 


rate  constants  of  1.62  x 10 


Hg.  The  cross  section  for  electron-impact  excitation  of  Hg  is 
limited  to  the  excitation  cross  section  for  the  Hg*  (3P2)  for  reasons 
discussed  in  Chapter  2.  Rockwood's  (Ref  58)  cross  sections  for  excitation 
to  various  Hg  states  are  the  result  of  generating  a self-consistent 
set  of  cross  sections  based  on  various  experimental  and  calculated 
transport  data.  Rockwood  (Ref  38)  uses  Borst  (Ref  59)  and  von  Engel 
(Ref  60)  for  his  experimental  data  on  excitation  of  Hg*  (3P2) 
Unfortunately,  none  of  the  three  authors'  cross  sections  covers  a 
large  range  of  energy.  The  theoretical  cross  section  of  McConnell  and 
Moiseiwitsch  (Ref  61)  is  used  in  this  calculation.  McConnell  and 
Moiseiwitsch  (Ref  61)  and  Rockwell's  (Ref  58)  calculations  are  in  good 
agreement.  3oth  find  a threshold  for  excitation  to  Hg*  (3P2)  at 
5.7  eV  and  a maximum  cross  section  of  3.2  x 10  ‘6  cm2  at  6.4  eV. 
McDonnell  and  Moiseiwitsch 's  cross  section  is  approximated  in  equation 
108. 


0 , 0 < E s 5.7,  e > 17  eV 

a = { 4.57145  x 10-15  e -2.60572  x lO-15  , 5.7  1 e 1 6.4  eV 

-0.5018S6S  x 10'15  e +5.1520755  x 10"16 , 6.4  i e i 17  eV  (10S) 

The  area  under  the  curve  generated  by  equation  108  is  within  1.5’o  of 
the  area  under  the  curve  of  McConnell  and  Moiseiwitsch  (Ref  61). 

Substitution  of  equation  108  into  equation  105  and  evaluation  over 
the  proper  limtis  yields  an  electron- impact  excitation  reaction-rate 
constant  of  8.45  x 10  cm3/sec  for  excitation  to  Hg*  (;P2) 
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Conclusion.  The  electron- impact  excitation  cross  sections  that 
are  calculated  in  this  section  using  a Maxwel i- Boltzmann  distribution 
for  electron  energies  may  be  incorrect.  The  Maxwell- Bolt zinann 
distribution  of  electron  energies  tends  to  skew  the  real  distribution 
of  electron  energies  toward  the  higher  energies.  Therefore,  the 
reaction-rate  constants  that  are  calculated  in  this  section  are 
probably  high. 


Ionization  of  Excited  States  by  Secondary-Electron  Impact 


Introduction.  The  ionization  process  by  secondary  electrons  is 


generally  in  the  form  of  reaction  109  where  X*  is  a metastable 
atom,  X+  is  the  ionized  atom,  and  e is  a secondary  electron. 


Q + x*  •*  X 


(109) 


This  process  is  a loss  term  for  metastables  and  a gain  term  for  secondary 
electrons  and  ions. 

The  ionisation  processes  that  are  included  are  the  ionization  of 
Ar*  , Xe*  , and  Hg*  . The  cross  sections  for  these  processes 
are  generated  using  equation  110.  Equation  110  is  an  empirical  expression 
for  ionization  cross  sections  of  excited  atoms  and  molecules  that  is 
developed  by  Vriens  (Refs  o2,  63). 
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The  cross  section,  a , is  a function  of  the  electron  energy  E , 

V 

in  eV , the  energy  necessary  to  ionite  the  excited  state  U , in 
eV  , and  nr  which  is  equal  to  0.S8  x 10  ' 0 cm2.  The  equation  is 
fairly  accurate  for  U < 3.5  eV  (Ref  62).  Grvzinski's  (Ref  64) 
expression  is  more  accurate  for  U > 10  eV  (Ref  62). 

The  cross  sections  that  are  generated  by  Vriens'  expression 
(Refs  62,  63)  are  approximated  by  straight  lines  as  in  equation  103. 

< These  expressions  are  substituted  into  equation  105  and  evaluated 
over  the  proper  limits  of  integration.  The  electrons  are  assumed 
to  be  given  by  a Maxwel 1- Boltzmann  distribution  with  4 eV  average 
energy. 

Ar*,  Xe*,  and  Hg*.  The  following  data  is  used  in  the  calculations 
for  the  ionization  of  Ar*  , Xe*  , and  ilg*  by  the  secondary 
electrons.  Assuming  the  Ar*  is  the  lowest  metastable  state,  the 
metastable  energy  is  11.543  eV.  Tlie  ionization  potential  for  ground 
state  Ar  is  15.75  eV  (Ref  65).  The  resultant  U is  4.211  eV.  The 
equation  describing  the  approximate  cross  section  for  Ar*  ionization 
is  given  in  equation  111. 

0 , 0 < e < 4.211  eV 

o = { 3. 68S2  x 10"16  e - 1.5531  x 10-"5,  4.211  < s < 6.1089  eV 
7 x 10“ 1,5  , s > 6.1089  eV  (111) 


Assuming  Xe*  is  in  the  lowest  metastable  state,  the  metastable  energy 
is  S.315  eV.  The  ionization  potential  for  ground  state  Xe  is 
12.129  eV  (Ref  65).  The  resultant  U is  3.314  eV.  The  equations 


describing  the  approximate  cross  section  for  Xe*  ionization  are 
given  in  equation  112. 

, 0 < e < 5.814  eV 

2.2761  x 10*15,  5.814  < e < 5.322  eV 
, e > 5.322  eV  (112) 


0 

a = { 5.9677  x 10_1 
9 x 10'16 


The  Hg*  is  considered  to  be  in  the  3P2  metastable  state  for  reasons 
discussed  in  Chapter  2.  The  excitation  energy  of  Hg*  (3Pn)  is 
5.454  eV,  and  the  ground  state  ionizaton  potential  of  Hg  is  11.434  eV 
(Ref  66).  The  resultant  u is  5.00  eV.  The  equations  describing 
the  approximate  cross  section  for  Hg*  ionisation  are  given  in 
equation  113. 


a 


0 , 0 < e < 5.0  eV 

= {2.48  x 10"16  - 1.24  x l(f15,  5.0  < e < 7.097  eV 

5.2  x 10'16  , e > 7.097  eV 


(113) 


Conclusions.  The  reaction-rate  constants  for  ionization  of 
Ar*  , Xe*  , and  Hg*  are  3.27  x 10~3  cm3/sec,  4.83  x 10"3  cm3/sec, 
and  1.94  x 10~3  cm3/sec,  respectively.  Again,  these  reaction-rate 
constants  are  probably  too  high,  because  the  Maxwell-Boltzmann 
distribution  is  used  to  describe  the  electron  energies. 

Bauer  and  Burtky  (Refs  t>7,  b8)  present  classical  methods  for 
calculating  cross  sections  for  both  excitation  and  ionization  by 
electron- impact . 


76 


Cl-Cl  Recombination  Reactions 


The  two  Cl-Cl  atom  recombination  reactions  included  in  this 
model  are  of  the  form  in  reaction  114. 


Cl  + Cl  + M -►  Cl2  + M CM  = Cl2,  Ar)  (114) 


The  reaction-rate  constants  for  this  model  are  the  result  of  experi- 
mental data  obtained  by  Clyne  and  Stedman  (Ref  69).  This  data  implies 
that  equation  115  gives  the  reaction-rate  constant  for  reaction  114 
with  Cl2  as  the  third  body.  The  temperature  T is  in  units  of  ° K. 


R 


10 


16  3 


.2  7 


c m bmole 


-2  -1 
sec 


(115) 


Subsituting  the  temperature  for  this  model,  533°  K,  into  equation  115 
results  in  a reaction-rate  constant  of  4.227  x 10' 5 cm6mole~2sec  * 
for  Cl2  as  a third  body.  This  reaction-rate  constant  converts  to 
1.165  x 10  32  cm°/sec 

The  data  from  Clyne  and  Stedman' s work  t.Ref  69)  implies  that 
equation  116  gives  the  reaction-rate  constant  for  reaction  114  with 
Ar  as  a third  body. 

-3 .0 

R = 1C)15'5  CsJgO  cm6mole-2sec-1  (116) 

Again,  temperature,  T , is  in  0 K.  Substitution  of  the  temperature 

for  the  model,  533°  K,  into  equation  116  results  in  a reaction-rate 

constant  of  7.099  x 101£*  cmsrnole  "sec  1 . This  reaction-rate 

-S'? 

constant  converts  to  1.957  x 10  cmc/sec  for  reaction  114  with  Ar 
as  the  third  body. 


The  value  found  for  reaction  with  Ar  as  the  third  body  by 

Clyne  and  Stedman  (Ref  69)  is  in  good  agreement  with  the  values  found 

by  Hutton  and  Wright  (Ref  70)  and  Bader  and  Ogrytlo  (Ref  71).  Clyne 

and  Stedman's  value  for  Ar  as  a third  body  in  reaction  114  is 

(4.4  ± 0.5)  x 101S  cmsmole"2sec"i  at  298°  K (Ref  69).  Hutton  and 

Wright's  value  for  the  same  reaction  at  300°  K is  (4.25  ± 0.8)  x 1015 
2 1 

* cm°mole'  sec-  (Ref  70).  Bader  and  Ogrvzlo's  value  for  this  reaction 

at  313°  K is  4 x 1015  cm6mole"2sec-1  (Ref  71). 

Radiative  Decay  of  ArCl*,  XeCl*,  and  HgCl* 

The  radiative  lifetimes  that  are  used  in  this  model  for  ArCl* 
and  XeCl*  are  based  on  experimental  data.  Gundel  et  al.  (Ref  56) 
estimate  the  radiative  lifetime  of  ArCl*  to  be  20  ns  at  higher  pressures. 
At  pressures  of  1 to  2 torr,  Gundel  et  al.  (Ref  56)  estimate  the  lifetime 
of  ArCl*  to  be  100  ns.  Hart  and  Searles  (Ref  72)  estimate  the  radia- 
tive lifetime  of  XeBr*  to  be  17.5  ns.  Dunning  and  Hay  (Ref  -3) 
calculate  the  radiative  lifetime  of  KrF*  to  be  6.5  ns.  These  radiative 
lifetimes  are  used  as  guidance  to  estimate  the  radiative  lifetimes  of 
ArCl*  and  XeCl*  to  be  20  ns.  These  radiative  lifetimes  lead  to 
reaction-rate  constants  of  5 x 107  sec-1  for  the  bound-free  radiative 
decay  of  ArCl*  (Ref  56)  and  XeCl*  to  their  products  of  Ar  and  Cl  - 
and  Xe  and  Cl  , respectively. 

The  radiative  lifetime  of  HgCl*  is  reported  to  be  40  ns  (Ref  S) . 

This  radiative  lifetime  for  HgCl*  leads  to  a reaction-rate  constant  of 
2.5  x 10'  sec-1  for  the  radiative  decay  of  HgCl*  . The  HgCl*  decays 
to  a bound  state  of  hgCl 
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Ionization  and  Excitation  by  Electron  Beam 

The  cross  sections  for  ionization  and  excitation  of  Ar  and  Xe 

by  electron-beam  electrons  are  provided  by  Hunter  (Ref  5) . These  cross 

sections  are  calculated  by  multiplying  2.4  times  the  Berger-Seltzer 

cross  sections  (Ref  74).  The  2.4  factor  is  a compensation  factor  for 

scattering  losses.  The  total  cross  section  is  then  partitioned  so  that 

approximately  75%  results  in  ionization  and  25%  results  in  excitation 

of  Ar  and  Xe  . For  the  beam  electron  energy  that  is  used  in 

Chapter  4,  the  cross  section  for  ionization  of  Ar  is  7.68  x 10“  cm2  ; 

the  cross  section  for  excitation  of  Ar  to  a metastable  state  is 
.I8 

2.11  x 10  cm2  ; the  cross  section  for  ionization  of  Xe  is 

1.90  x 10' 1 cm2  ; and  the  cross  section  for  excitation  of  Xe  to  a 
metastable  state  is  6.33  x 10"  cm2  . Other  work  on  ionization 


sections  for  Ar  is  provided  in  References  75  and  76,  and  does  agree 
fairly  well  with  Berger  and  Seltzer  (Ref  74) . 


HgCl*  Quenching 

The  only  reaction  that  is  included  in  this  model  for  the  quenching 
of  HgCl*  is  reaction  32. 


HgCl*  + Cl2  - HgCl  + Cl2 


The  reaction-rate  constant  for  this  reaction  is  similar  to  that  for  XeBr* 

quenching  by  Br2  and  KrF*  quenching  by  F2  . The  reaction-rate 

.10  , 

constant  for  XeBr*  quenching  by  Bt2  is  (S  t 2)  x 10  cnWsec 
(Ref  "2).  The  reaction-rate  constant  for  KrF*  quenching  by  F2  is 
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.9  , 

1 x 10  cnr/sec  CRef  77).  The  reaction-rate  constant  for  HgCl*  by 

-10  , 

CI2  is  reported  to  be  5.62  x 10  cnr/sec  (Ref  8).  This  reported 
value  is  the  quenching  reaction-rate  constant  used  in  this  model. 

The  difference  between  quenching  of  both  XeBr*  by  Br2  and  KrF* 
by  F2  and  HgCl*  by  CI2  is  that  a bound  state  exists  for  HgCl  , 
but  does  not  exist  for  XeBr  and  KrF  , All  HgCl*  appears  to  be 
quenched  to  a bound  state  of  HgCl  (Ref  8) . 

Summary 

The  reaction-rate  constants  that  are  developed  in  this  section  are 
based  on  both  theoretical  and  experimental  data.  Table  2 shows  the 
reactions  and  reaction-rate  constants  that  are  used  in  the  solution  of 
this  model  that  is  presented  in  Chapter  4. 

Table  II 

Reactions  § Reaction-Rate  Constants 


Reaction 

Equation 

Number 

Reference 

Reaction-Rate 

Constant 

e + 

Ax*  0 

Ar*  + Ar 

Text 

17 

3.9  x 10  3 cm3/sec 

e + 

Xe£  -> 

Xe*  + Xe 

Text 

18 

8.16  x 10"3  cm3/sec 

e + 

ci2  - 

Cl"  + Cl 

Text 

23 

3.7  x 10  " cm3/sec 

Ar+ 

+ 2Ar 

-*  Ar2  + Ar 

Text 

31 

-31 

2.5  x 10  cmb/sec 

Xe+ 

+ Xe 

+ Ar  -*•  Xe£  + 

Ar 

77b 

27 

4.96  x 10"'1  cm6/sec 

+ 

Xe 

+ Xe 

+ Xe  -►  Xe2  + 

Xe 

Text 

31 

3.6  x 10  cm6/sec 

Ar2 

+ Xe 

Xe+  + 2Ar 

S3 

36 

-10  , 

6.65  x 10  cnr/sec 

Xe2 

+ Hg 

-v  Hg+  + 2Xe 

83 

36 

-10  , 

3.11  x 10  cnr/sec 

Ar* 

+ Cl* 

- ArC 1 * 

97a 

33 

2.13  x 10"°  cm3/sec 

Xe" 

+ Cl" 

- XeCl* 

97a 

33 

1,75  x 10  ° cm3/sec 
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Table  II  (Cont’d) 


Reaction 


Hg+ 

+ Cl 

' HgCl* 

Hg* 

+ Cl 

2 HgCl*  + Cl 

Xe* 

+ Cl 

2 -*■  XeCl*  + Cl 

Ar* 

+ Cl 

2 -*■  ArCl*  + Cl 

e + 

Ar  -» 

• Ar*  + e 

e + 

Xe  -> 

■ Xe*  + e 

e * 

Hg  -> 

• Hg*  + e 

e + 

Ar* 

-*■  Ar+  + 2e 

e + 

Xe* 

-*■  Xe+  + 2e 

e + 

Hg* 

-*•  Hg+  + 2e 

Cl  + 

Cl 

+ Cl2  2C12 

Cl  ♦ 

• Cl 

+ Ar  Cl2  + Ar 

ArCl 

* 

Ar  + Cl  + hv 

XeCl 

.*  -*■ 

Xe  + Cl  + hv 

HgCl 

.* 

HgCl  + hv 

P . + P 

e * At  -*■  Ar  + e + e 

e^  + Ar  -*■  Ar*  + e 

e + Xe  -*■  Xe  + e + e 

P P 

e + Xe  -*•  Xe*  + e 

HgCl*  + Cl2  - HgCl  + Cl2 


Equation 

Number 

Reference 

97a 

33 

Text 

5 

Text 

52 

Text 

54,55 

105,106 

57 

105,106 

57 

105,108 

61 

105,111 

62 

105,112 

62 

105,113 

62 

115 

69 

116 

69 

Text 

56 

Text 

Text 

Text 

8 

Text 

5 

Text 

5 

Text 

5 

Text 

5 

Text 

8 

Reaction-Rate 

Constant 

1.685  x 10~6  cm3/sec 

-10  5 

4.18  x 10  cnr/sec 

1.04  x 10'  cm3/sec 

-10  , 

3.16  x 10  cnr/sec 

-10  , 

1.62  x 10  cm3/sec 

6.92  x 10'10  cm3/sec 

8.45  x 10  " cm3/sec 
-8  , 

3.27  x 10  cnr/sec 
-8  , 

4.85  x 10  cnr/sec 
-8  , 

1,94  x 10  cnr/sec 
-32 

1.165  x 10  cnr/sec 

1.957  x 10  cm°/sec 

, -1 
5 x 107  sec 

5 x 107  sec'1 

2.5  x 107  sec" 


- 1 0 , 

5,62  x 10  cnr/sec 


The  reaction-rate  constants  for  the  excitation  of  Ar  , Xe  , and 
Hg  and  the  ionisation  of  Ar*  , Xe*  , and  Hg*  by*  secondary  electrons 
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may  need  to  be  adjusted  in  Chapter  4,  The  assumption  that  the  secondary 
electrons  are  in  a Maxwell -Bo lttmann  distribution  with  an  average 
energy  of  4 eV  is  probably  incorrect. 
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IV,  Solution  and  Analysis 


Overview 

In  this  chapter,  the  rate  equations  for  the  various  species  that 
are  presented  in  Chapter  2 are  used  to  determine  solutions  for  two 
cases  of  this  model.  The  first  case  that  is  investigated  involves  the 
reaction  kinetics  of  a HgCl  laser  that  is  generated  by  an  electron 
beam  only.  The  second  case  involves  the  reaction  kinetics  of  a HgCl 
laser  that  is  initiated  by  an  electron  beam,  and  the  lasing  plasma  is 
subjected  to  an  applied  d.c,  electric  field. 

A computer  program  supplied  by  A.  M.  Hunter  (Ref  5j  is  used  to 
analyze  both  cases  for  various  electron  beam  currents  and  reaction-rate 
constant  sensitivities.  This  program  uses  a numerical  integrator  to  solve 
the  rate  equations,  The  integrator  solves  for  the  change  in  number  densities 
of  the  various  species  by  stepping  forward  in  time  by  a small  time  interval 
and  adding  or  subtracting  the  increase  or  decrease  in  number  densities 
to  the  appropriate  number  densities  for  each  species  at  the  beginning  of 
the  small  time  interval. 

Initial  computer  program  runs  did  not  indicate  a stable  plasma  could 
be  maintained  for  the  discharge-sustained  HgCl  laser  simulation.  This 
problem  is  caused  by  the  inaccurate  assumption  that  the  secondary  electrons 
can  be  described  by  a .Maxwell -Boltzmann  energy  distribution.  This  problem 
is  solved  by  adjusting  the  reaction-rate  constants  for  excitation  of  Ar 
and  Xe  by  secondary  electron  impact  and  for  ionization  of  Ar*  and 
Xe*  by  secondary  electron  impact,  By  adjusting  these  reaction-rate 
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constants  to  the  constants  used  for  the  same  reactions  in  the  KrF  laser, 
a stable  plasma  is  achieved  for  the  discharge-sustained  HgCl  laser 
simulation.  The  constants  that  are  used  for  the  solution  in  this  chapter 

-12  i 

are  6.62  x 10  cmJ/sec  for  excitation  of  Ar  by  secondary  electrons, 

7.46  x 10  cm3/sec  for  ionization  of  Ar*  by  secondary  electrons, 

- 1 1 -J 

7.75  x 10  cm-Vsec  for  excitation  of  Xe  by  secondary  electrons, 
and  1.16  x 10  7 cm3/sec  for  ionization  of  Xe*  by  secondary  electrons. 

I 

These  constants  are  provided  by  A.  M,  Hunter  (Ref  51  and  eliminate  the 
problems  for  these  reactions  that  are  caused  by  assuming  a Maxwell- 
Boltzmann  distribution  for  the  secondary  electrons.  The  Maxwell -Boltzmann 
distribution  contains  more  higher  energy  secondary  electrons  than  the 
actual  secondary  electron  energy  distribution  (Ref  51 . 

This  chapter  is  developed  in  the  following  manner.  First,  a 
determination  is  made  to  find  which  species  number  densities  vary  slowly 
with  time  and  can  be  considered  constant.  Second,  an  analytical 
solution  of  this  model  with  no  applied  electric  field  is  presented  for 
the  electron  density  and  is  compared  to  the  computer  solution.  Third, 
an  analytical  solution  of  this  model  with  an  applied  electric  field  is 
presented  for  the  electron  density  and  is  compared  to  the  computer 
solution.  The  primary  emphasis  is  placed  on  the  electron  number  density, 
because  the  electron  number  density  is  found  to  determine  the  stability 
of  the  e-beam  sustained  KrF  laser  (Ref  5)  and  is  assumed  to  determine 
the  stability  of  the  HgCl  laser  in  a similar  manner. 


Slow  Varying  Species 

Before  an  analytical  solution  is  found  for  this  model,  the  species 
number  densities  that  vary  slowly  are  determined.  A species  number 
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density  that  is  found  to  vary  slowly  can  be  assumed  to  be  constant  in 
the  analytical  solution. 


A species  number  density  is  considered  to  vary  slowly  if  the  initial 
lifetime  for  an  atom  or  molecule  of  a species  is  greater  than  10-t> 
seconds.  The  average  initial  lifetime  for  an  atom  or  molecule  is 
determined  by  dividing  the  initial  n^  by  the  initial  number  density 
of  the  species.  The  resultant  number  is  approximately  equal  to  the 
inverse  of  initial  lifetime.  This  relationship  is  shown  in  equation  117. 


(117) 


The  only  species  that  exist  at  time  zero  are  Ar  , Xe  , Hg  , Cl2  , 

p 

and  e . All  other  species  number  densities  are  initially  zero. 
Therefore,  the  only  species  with  an  initial  lifetime  are  Ar  , Xe  , 

Hg  , Cl2  , and  electrons. 


The 

Ar 

, Xe 

. Hg 

are  found 

in 

the  following 

equations 

35, 

36,  37, 

44,  ; 

n.  for  Ar 

l 

, Xe  , 

Hg 

4.531  x 1C 

>13 

cm" 3 for 

Xe  , 

cm-3  for 

Cl2 

. The 

n. 

l 

is  the  number 

density 

for 

lifetimes 

for 

Ar  , 

Xe 

-3 


iia  ,~m-3 


10  seconds  for  electron  beam  currents  to  15  amps/cnr  . The  initial 
lifetime  for  electrons  is  found  to  be  less  than  10'6  seconds  for  all 
electron  beam  currents  greater  than  1 amp 'em2  . Therefore,  the  number 
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densities  for  Ar  , Xe  , Hg  , and  CI2  are  considered  constant 
in  the  analytical  solution.  The  number  density  of  electrons  is  not 
considered  constant  in  the  analytical  solution. 

Solution  for  HgCl  Laser  Without  Applied  Electric  Field 

Introduction  and  assumptions.  The  analytical  solution  for  the 
electron  number  density  is  based  on  two  major  assumptions . One 
smaller  assumption  is  made  in  the  development  of  this  solution.  This 
solution  is  done  for  a 5 amp/cm2  electron  beam  current  and  an  initial 
gas  mixture  of  3.726  x 1019  Ar  atoms/cm3  , 4.531  x 1019  Xe  atoms/cm'  , 
1.745  x 1013  Hg  atoms/cm3  , and  9,062  x 1016  CI2  molecules/cm3 
The  reaction-rate  constants  that  are  used  in  this  solution  are  those 
reaction-rate  constants  that  are  presented  to  this  point  in  this  paper. 

The  primary  assumption  that  is  made  for  the  analytical  solution  of 
the  HgCl  laser  without  an  applied  electric  field  is  that  the  secondary 
electrons  are  in  thermal  equilibrium  with  the  surrounding  gas.  This 
simplifying  assumption  is  not  completely  true,  because  it  does  not 
include  the  effect  of  thermalizing  the  secondary  electrons  that  start 
at  a higher  than  thermal  energy.  A true  distribution  of  these  secondary 
electrons  can  be  found  by  solving  the  Boltzmann  equation,  but  is  beyond 
the  scope  of  this  paper.  This  assumption  of  thermal  electrons  changes 
the  dissociative-recombination  reaction-rate  constants  for  Art  (reaction 
10)  and  Xet  (reaction  21)  to  5.99  x 10"'  cm3/sec  and  1.S8  x 10*bcm3/sec 
respectively  (Refs  17,  18).  These  increases  are  due  to  the  decrease  in 
secondary  electron  temperature. 
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A second  assumption  that  is  made  in  the  analytical  solution  of 
the  HgCl  laser  without  an  applied  electric  field  is  that  the 
secondary’  electrons  are  not  energetic  enough  to  excite  Hg  , Ar  , or 
Xe  to  their  metastable  states  and  are  not  energetic  enough  to  ionite 
the  excited  states  of  Hg*  , Ar*  , and  Xe*  . Therefore,  the 
reaction-rate  constants  for  reactions  15  through  20  are  set  equal  to  zero. 
This  assumption  is  not  totally  correct,  because  the  effect  of  thermalizing 
the  secondary  electrons  that  start  at  higher  than  thermal  energies  is 
not  included  in  this  analysis. 

The  third  assumption  that  is  made  in  this  solution  is  that  an 
equilibrium  electron  number  density  exists  for  at  least  1 x 10'6 
seconds  and  that  equilibrium  number  densities  exist  for  all  other  species. 
This  assumption  is  almost  accurate,  because  these  number  densities  do 
in  fact  vary  slowly  after  they  are  initially  established  for  all  species. 
Also,  the  equilibrium  electron  number  density  is  assumed  to  be  much  less 
than  5.0  x 1015  electrons/cm3  . 

Analytical  solution.  The  analytical  solution  is  a result  of  the 

solution  of  equations  35,  38,  59,  40,  and  41.  The  equilibrium  number 

density  of  electrons  is  found  by  setting  ng  equal  to  zero  and  solving 

for  n . The  n is  found  to  be  dependent  on  nv  + and  n + 
e e r Xea  Ar2 

By  setting  equation  41  equal  to  zero,  n,  + is  found  to  be  equation  1 1 S , 

if  n,  + and  n are  in  cm'3 
Ar  e 


n 


Art 


C3.471  x lO8)  n.  + , 

— cm' 6 Cl  IS) 

C3.013  x lO3)  + (5.99  x 10"  ) ng 


87 


AIR  FORCE  INST  OF  TECH  WRIGHT -PATTERSON  AFB  OHIO  SCH— ETC  F/G  20/5 
REACTION-R  CONSTANTS  AND  SIMPLE  RATE  THEORIES  APPLIED  TO  THE  ME— ETC (U) 
MAR  78  R P SUMMERS 


AF  TT/6NE/PH/78M— 9 


NL 


If  the  equilibrium  electron  number  density  is  assumed  to  be  less  than 
5.0  x 1015  cm"3  , equation  118  reduces  to  equation  119, 


n,  + 
Ar2 


0.1152  n.  + cm 
Ar 


(119) 


The  solution  for  the  equilibrium  Ar  number  density  is  given  in 

-3 


equation  120.  The  nCj-  is  in  cm 


Substitution  of  equation  120 
into  equation  119  results  in  equation  119a. 


(8.942  x 1021) 


n.  + 
Ar 


-6> 


cm 


.3 


(120) 


(3 . 471  x 10a)  + (2,134  x 10  ) nCJ- 


(1.03  x 1021) 


n,  + 
Ar2 


cm 


- 3 


(119a) 


(3.999  x 107)  + (2.458  x 10'7)  nci- 


The  equilibrium  number  density  of  Xe2  is  found  by  solving  equation 
38  for  the  equilibrium  Xe+  number  density  and  substituting  the  Xe+ 
solution  into  the  Xe2  equation,  Equation  119a  is  substituted  into 
the  Xe2  solution  to  give  equation  121. 


(9.113  x 107) 


Xe2  (5.427  x 108)  + (1.88  x 10"6)  n 


(121) 


(3.211  x 1030)  ♦ (6.61  x 1014)  ncl- 


(3.644  x 1015)  + 92.38  ncl.  + (4.3  x 10"13)  ncl.2 
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Equations  121  and  119a  are  substituted  into  equation  35  to  give  the 
quadratic  equation  122  for  ng  . 


(1.16  x 109) 

[(6.293  x 102)  - ] n2 

C3.999  x 107)  ♦ (2.458  x 10-7)  na_ 


(3.348  x 1023) 

[(2.037  x 1017)  - 

(3.999  x 107)  + (2.458  x 10"  ) n r 


(3.211  x 103°)  ♦ (6.61  x 1014)  n 

(1.71  x 102)  { — — > ] 

(3.644  x 1015)  ♦ 92.38  n ♦ (4.3  x 10'13)  nc|_  e 


(6.313  x 1030) 

0 


(122) 


The  solution  of  the  quadratic  equation  122  yields  an  equilibrium  ng  of 
5.683  x 1013  cm'3  for  an  ncl-  of  6.488  x 1013  cm'3  . This  solution 
agrees  well  with  the  computer  results  that  are  shown  in  Table  3. 

Equation  122  also  results  in  a decrease  in  equilibrium  ng  as 
nCj_  is  increased  and  an  increase  in  equilibrium  ng  as  n^-  is 
decreased.  Examination  of  equation  48  shows  that  the  only  source  of 
n_._  is  the  dissociative  attachment  of  CI2  by  secondary  electrons. 
Therefore,  the  equilibrium  ng  is  sensitive  to  the  reaction-rate  constant 
for  dissociative  attachment  of  CI2  as  shown  in  Table  3. 
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Table  III 

Equilibrium  ng  as  a Function  of  n^. 
Beam  Current  is  5 amp/ cm2. 


Dissociative  Attachment 
Reaction-Rate  Constant 
for  Cl2  (cm3/sec) 


Equilibrium  n^. 
Population  (cm-3) 


Equilibrium  ng 
Population  ('em-3) 


1.0 

X 

CD 

1 

o 

H 

7.717 

X 

1013 

1.439 

X 

1013 

3.7 

X 

10'9 

6.488 

X 

1013 

3.809 

X 

1013 

1.0 

X 

10-10 

4.889 

X 

1012 

5.487 

X 

1014 

As  is  expected,  the  equilibrium  ng  increases  as  the  initiating 
beam  current  is  increased.  This  is  shown  in  Table  4 where  only  the 
electron  beam  current  is  varied. 


Table  IV 

Equilibrium  ng  as  a Function  of  Electron  Beam  Current 
Computer  Solution  at  1 x 10"6  Seconds 


Electron  Beam  Current 
(amp/ cm2) 

1 

2 

S 

10 


Equilibrium  ng 
Population  (cm"3) 

7.032  x 1012 

1.394  x 1013 

3.809  x 1013 

7.891  x 1013 


Conclusion.  The  equilibrium  number  density  of  electrons  is  approxi- 
mately proportional  to  the  incident  electron  beam  current.  This 
equilibrium  ne  is  also  fairly  sensitive  to  the  dissociative  attachment 
rate  for  Cl2  as  is  shown  in  Table  3. 


Solution  for  the  HgCl  Laser  with  Applied  Electric  Field 


Introduction  and  assumptions.  The  analytical  solution  for  the 
electron  number  density  for  the  HgCl  laser  with  an  applied  electric 
field  is  done  for  a beam  current  of  1 amp/cm-  and  the  same  number 
densities  that  are  used  in  the  solution  for  the  HgCl  laser  without 
an  applied  electric  field.  Again,  the  reaction-rate  constants  that 
are  used  in  this  solution  are  those  reaction-rate  constants  that  are 
presented  to  this  point  in  this  paper.  The  assumptions  for  this 
solution  are  different  from  the  assumptions  that  are  used  in  the 
solution  for  the  HgCl  laser  without  3n  applied  electric  field. 

The  first  assumption  is  that  the  electric  field  maintains  the 
secondary  electrons  at  an  average  energy  of  4 eV.  This  number  is 
similar  to  the  energy  of  secondary  electrons  in  the  KrF  laser  (Ref  5). 
This  assumption  changes  the  dissociative  recombination  reaction-rate 
constants  for  Art  (reaction  10)  and  Xej  (reaction  21)  to  the 
reaction-rate  constants  for  these  reactions  that  are  presented  in 
Table  2.  This  assumption  also  means  that  enough  energy  now  exists  in 
some  secondary  electrons  to  cause  excitation  to  metastable  states  of 
Hg  , Ar  , and  Xe  . Also,  some  secondary  electrons  can  ionize  Hg* 
Ar*  , and  Xe*  . The  reaction-rate  constants  that  are  used  for  Hg 
excitation  and  Hg*  ionization  are  the  reaction-rate  constants  that 
are  presented  in  Table  2.  The  reaction-rate  constants  that  are  used 
for  Ar  and  Xe  excitation  and  Ar*  and  Xe*  ionization  by  secondary 
electrons  are  the  reaction-rate  constants  that  are  presented  in  the 
overview  of  this  chapter. 
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The  second  assumption  is  that  an  equilibrium  number  density  for 
all  species  does  exist  at  some  time.  The  electron  equilibrium  number 
density  is  assumed  to  be  within  an  order-of -magnitude  of  the  electron 
equilibrium  number  density  that  is  found  for  the  HgCl  laser  without 
an  applied  electric  field.  This  assumption  is  probably  a reasonable 
one,  because  the  same  species  number  densities  are  being  used,  3nd 
Ar  , Hg  , CI2  , and  Xe  still  all  very  slowly. 

Analytical  solution.  This  solution  is  based  on  equation  35,  the 
electron  rate  equation.  The  first  step  is  to  investigate  the  loss  terms 
of  equation  35  to  see  if  a dominant  loss  term  is  found.  These  terms 
are  listed  in  equation  123. 


- R«  n n„,  - R,rt  n n + - R.,,  n nv  + 
9a  e Cl2  10  e Ar2  21  e Xe 


(123) 


The  ng  factor  is  a common  factor  to  all  terms  in  equation  123  and  can 
be  ignored  for  the  present  argument.  By  approximating  R.^  as  2R^0 
and  eliminating  the  common  n^  factor,  equation  123  is  approximately 
equation  124. 


R9a  nCl2 


"Xe^ 


(124) 


Because  the  entire  plasma  maintains  overall  neutral  state  and  charge 
is  conserved,  the  number  of  negative  charges  is  equal  to  the  number  of 
positive  charges.  This  relationship  is  shown  in  equation  125. 


r 


n ♦ n„ , . 
e Cl 


n.  + ♦ nv  + ♦ n„  ♦ ♦ n,  ♦ + nv  ♦ 

Ar  Xe  Hg  Ar2  Xe2 


(.125) 


The  number  of  Cl2  varies  slowlv,  Therefore,  nri  is  much  greater 

*-  * 2 

than  n„,_  . The  sum  of  n,,,_  and  n is  much  less  than  n...  , 

Li  ci  e c i2 

because  ng  is  probably  about  1013  cm-3  , and  n^.^  is  almost 
1017  cm"3  . Equation  126  shows  the  relationship  between  n^  and 
both  nri_  and  n . 

LI  © 


n_,  >>  n^  + n_._ 

Cl?  e Cl 


Solving  equation  125  for  the  sum  of  n,  ♦ and  nv  ♦ and  using 

A12  a ©2 

equation  126  results  in  equation  127. 


(126) 


n,  «■ 
\r2 


nCl2  >>s*  ne  * nCl" 


nXe*  * nHg* 


n,  ♦ > 

Ar 


♦ nv  ♦ 
Xe? 


(127) 


The  R„  reaction-rate  constant  is  onlv  an  order-of-magnitude  less  than 
9a 

R,-  , and  n_.  is  several  orders-of -magnitude  greater  than 

10  L 1 2 

nAr2 

(— * — * nY  + ) . Therefore,  equation  128  is  reasonable  and 

Z \©2 

R„  n n„,  is  the  dominant  electron  loss  mechanism. 

9a  e Cl2 


nAr? 

R9a  ne  nCl2  ~ 2R10  * "xe? 
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(128) 


. r 


Equation  35  is  now  approximated  by  equation  35a. 


n » V .a,  n ,n,  ♦ V .crv  n ,nv  ♦ R,,n  n„  * 

e eb  Ar  eb  Ar  eb  Xe  eb  Xe  16  e Hg* 


♦ R,0n  n.  . ♦ R^rtn  nv  * - Rn  n n 

18  e Ar*  n » q->  ^ 


(35a) 


20  e Xe* 


9a  e Cl; 


The  equilibrium  n&  solution  is  found  by  setting  equation  35a 


equal  to  ;ero  and  solving  for  ng  as  is  shown  in  equation  129 


V ,o.  n . n,  ♦ 
eb  Ar  eb  Ar 

ne  ~ R.  n 


- ft",  *11. 


VebaXenebnXe 


9a“Cl2  ' '■,'16‘‘Hg*  * R18nAr*  * R',r'n 


20  Xe* 


) 


(129) 


The  denominator  of  equation  129  is  the  determining  factor  for  the 


equilibrium  n^  . By  assuming  that  R^  n^  is  the  dominant 


factor  in  the  denominator  and  the  other  terms  in  the  denominator  can 


be  ignored,  the  equilibrium  n^  for  an  electron  beam  current  of 


1 amp/cnr  is  6.937  x 1012  electrons  cm"3  . This  number  is  within  a 


factor  of  1.5  of  the  number  that  is  found  by  using  the  computer  program 
at  1 amp/ cm2 


Analvsis  of  equation  129  shows  that  as  R„  is  varied  the  n 
1 9a  e 


grows  larger  or  smaller.  If  R^a  is  smaller  than  the  value  that  is 


reported  by  Sides  et  al.  (Ref  23),  the  equilibrium  n&  is  greater. 


If  Rga  is  larger  than  the  value  that  is  reported  by  Sides  et  al. 


(Ref  23),  the  equilibrium  n^  is  smaller,  If  R is  so  small  that 


9a 


the  denominator  is  almost  :ero,  no  equilibrium  n&  exists.  Therefore, 


the  dissociative-attachment  reaction-rate  constant  for  Cl?  at  electron 


!■ 
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energies  about  4 eV  needs  to  be  determined.  The  effects  of  various 


Rga  values  is  shown  in  Table  5 for  the  HgCl  laser  with  an  applied 
electric  field. 

Table  V 

Computer  Solution  for  HgCl  Laser  and  R^a  Varied 
1 amp/ cm-  Applied  Electric  Field 

■j 

3 Discharge  Transition 

R^a  (cm3/sec)  Equilibrium  n^  (cm*  ) fnsec) 


1.0 

X 

* 0 

10  • ' 

S.672 

X 

1013 

32.3 

3.7 

X 

© 

1 

10 

1.172 

X 

10*  3 

224 

1.0 

X 

10'9 

4.775 

X 

10*’ 

> 1000 

Tables  shows  that,  as  the  equilibrium  n increases,  the  plasma 
stability  is  more  difficult  to  maintain.  The  discharge  transition 
column  in  Table  5 indicates  the  time  at  which  the  plasma  is  definitely 
unstable  and  lasing  action  is  destroyed. 

A similar  effect  to  that  of  the  HgCl  laser  without  an  applied 
field  is  found  for  the  HgCl  laser  with  an  applied  electric  field 
as  the  electron  beam  current  is  increased.  The  effect  is  shown  in 
Table  6. 

Table  VI 

Computer  Solution  for  HgCl  Laser  with 
1 amp/cm3  Applied  Electric  Field  and  Electron  Beam  Cui'rent  Varied 

Electron  Seam  Current  (amp/ cm2)  Equilibrium  n?  (cm'3)  Discharge  Transition 

(nsec) 


l 

1.172 

X 

10* 3 

224 

2.180 

X 

10* 3 

93.7 

5 

6.122 

X 

10* 3 

41.4 
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Conclusion.  Again,  the  dissociative-attachment  reaction-rate 
constant  is  a determining  factor  for  the  equilibrium  n^  . The 
applied  electric  field  raises  the  equilibrium  ng  3nd  causes  the 
plasma  to  become  unstable  sooner  by  the  addition  of  electrons  that 
are  created  by  ionization  of  metastable  states. 

Summary 

The  dissociative-attachment  reaction-rate  constant  for  Cl2 
is  an  important  determining  factor  in  both  the  equilibrium  n^  and 
the  plasma  stability  time.  The  electron  population  appears  to  be 
an  important  factor  in  the  overall  plasma  stability. 
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V.  Conclusions  and  Recommendations 

Several  factors  concerning  the  model  that  is  constructed  and 
analyzed  in  this  paper  are  interesting  for  different  reasons.  The 
reaction-rate  constants  that  are  developed  in  Chapter  5 are  examples 
of  the  excellent  theoretical  and  experimental  predictions  that  can 
be  made  by  using  simple  theories  and  experimental  data.  Theoretical 
reaction-r3te  constants  are  predicted  for  ion-neutral  association 
reactions,  charge  transfer  reactions,  ionic  recombination,  and  har- 
pooning reactions.  The  theoretical  predictions  give  reasonable  results. 

The  reaction-rate  constants  that  are  presented  using  experimental 
data  as  a basis  are  sometimes  incomplete  or  in  need  of  further 
research.  For  example,  the  ArCl*  formation  by  Ar*  being  quenched 
by  Cl2  only  includes  501  of  the  total  quenching  of  Ar*  by  Cl2 
(Ref  56).  The  radiative  lifetimes  for  ArCl*  and  XeCl*  are  only 
estimates  that  are  based  on  similar  radiative  lifetimes.  The  ArCl* 
and  XeCl*  radiative  lifetimes  are  probably  within  a factor  of  three 
of  the  actual  radiative  lifetimes. 

Further  research  needs  to  be  conducted  into  the  dissociative- 
attachment  reaction-rate  constant  for  Cl2  at  higher  electron  energies. 
The  equilibrium  electron  number  density  and  HgCl  laser  stability 
are  sensitive  to  this  reaction-rate  constant  3t  all  electron  energies. 
The  majority  of  research  for  this  reaction  is  done  for  thermal 
electrons.  The  reaction- rate  constant  that  is  used  in  this  model  is 
for  thermal  electrons  (Ref  23)  and  is  probably  inaccurate  for  the 
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HgCl  laser  analysis  that  is  done  for  an  applied  electric  field. 

The  assumption  of  a secondary  electron  energy  distribution  that 
is  Maxwellian  is  incorrect.  This  assumption  is  made  to  simplify 
calculations.  Many  of  these  calculations  are  adjusted  in  Chapter  4. 

The  adjustments  that  are  made  in  Chapter  4 show  that  the  assumption 
of  a Maxwel 1-Boltzmann  distribution  for  secondary  electron  energy 
forces  to  the  reaction-rate  conztan^s  that  are  calculated  in  Chapter  3 
for  excitation  to  be  high.  The  reaction-rate  constants  that  are 
calculated  in  Chapter  3 for  ionization  of  metastable  states  by 
secondary  electrons  are  low.  Corrections  are  only  made  for  excitation 
of  Ar  and  Xe  and  ionization  of  Ar*  and  Xe*  by  secondary 
electrons.  The  excitation  of  Hg  and  ionization  of  Hg*  by  secondary 
electrons  needs  to  be  investigated  further.  A Boltzmann  equation 
solution  for  these  reactions  is  one  possible  approach.  This  approach 
possibly  can  be  used  to  investigate  primary  electron  interaction  with 
Hg  . 

One  aspect  of  high-energy  lasers  that  is  not  investigated  in 
this  paper  is  the  initial  distribution  of  vibrational  states  of  the 
primary  electron  energy  absorbing  gas.  A possible  dependence  exists 
between  the  gas  temperature,  the  initial  vibrational  state  distribution, 
and  several  reaction-rate  constants. 

A final  recommendation  for  future  research  concerns  the  effect 
of  pressure  on  thermolecular  charge  transfer  reactions.  These 
reactions  are  of  the  general  form  that  is  shown  in  reaction  130 
where  X and  Y are  different  atoms. 
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X,  ♦ Y,  + X - Y,  ♦ 5X 


(1301 


Recent  research  in  this  area  is  presented  in  References  "3,  "9,  and  30. 
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APPENDIX  A 


Other  Possible  Reactions  in  the  HgCl  Laser 


Xe* 

+ Xe  ♦ 

Ar  -*■  ArXe* 

+ Xe 

Ar* 

+ Ar  + 

Xe  -»  ArXe* 

+ Xe 

ArXe* 

+ Xe 

-*•  Xej  + Ar 

ArXe+ 

+ Ar 

-*>  Arj  + Xe 

Xej 

+ Hg  + 

Xe  -*■  Hg*  + 

3Xe 

Xe+ 

♦ Clj  - XeCl*  + Cl 

Ar+ 

+ Clj  ->•  ArCl*  ♦ Cl 

Hg+ 

+ Clj  -►  HgCl*  ♦ Cl 

O 

hv(5576  A)  + 

Ar**  -*•  Ar* 

+ e 

(Ref  81) 

hv(55 

76  A)  + 

Xe**  -*>  Xe* 

♦ e 

(Ref  81) 

e + 

Cl2  - 

Clj 

(Ref  11) 

Hg* 

+ Hg  + 

Hg  -*>  Hg$  + 

Hg 
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Photodetachment  of  Cl 

The  photodetachment  of  C]  is  shown  below. 

hv  ♦ Cl  •»  Cl  ♦ e 

This  reaction  is  not  included  in  this  model  because  the  cross  section 
is  zero  at  the  557.6  nm  wavelength  of  the  HgCl  laser  (Ref  32). 

This  zero  cross  section  for  the  photodetachment  of  electrons  from 
Cl*  indicates  the  HgCl  photons  are  not  at  a higher  energy  than 
the  threshold  energy  that  is  located  3t  a wavelength  of  approximately 
545.0  nm  for  this  reaction  (Ref  82). 
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